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FOREWORD 


This RECORD will be of primary interest to engineers and researchers 
who are involved with the underlying problems associated with roadway 
failures. A fewof the many ramifications of the frost problem are covered 
in the first 4 papers. Akili describes the stress-strain behavior of a 
clayey silt anda clay as evidenced by a brittle fracture or plastic deforma- 
tion. His study leads to a conclusion that the liquid water in the thin-film 
form controls the reaction. Wang and Roderick enhance the basic knowl- 
edge of frost-action reinforcing concepts that frost heave and ice segrega- 
tion are dependent on the initial (compacting) moisture, soil texture, and 
dry density. Laba studied the reactions of a sand-ice system subjected to 
temperature change below freezing. Jumikis contributes influence values 
for the 3-dimensional disturbance distribution for a heated rectangular 
structure on permafrost. Allof these expand our knowledge of frost action 
in soils. 

The remaining 3 papers are directed toward moisture problems in soils 
and bases in roadway structures. Haliburton, reporting on a 6-year study, 
gives recommendations regarding staged construction, geometric profile, 
shoulders, drainage, subbase, base, surfacing, and subgrade testing in 
relation to prevention of swell through proper consideration of the mois- 
ture problem. Healy and Long describe the construction and laboratory 
and field testing of a prefabricated drain concept using synthetic cloth and 
a channelized drainage core. McLean and Krizek describe the application 
of a finite-element method solution of a steady-state gravity flow to the 
design of subsurface drains for a depressed highway profile. 


CONTENTS 


STRESS-STRAIN BEHAVIOR OF FROZEN FINE-GRAINED SOILS 
Minden aie ss 5 ee oe evsiei he Gales Gi Blt or eae ee Aen Geo em dee Se 


FROST BEHAVIOR OF COMPACTED SOILS 
Mian=Chane Wang and Gilbert Li. Roderick. «66.006 ssw eww see ame we 
VISCOELASTIC PROPERTIES OF A LATERALLY CONFINED 


SAND-ICE SYSTEM SUBJECTED TO TEMPERATURE 
INCREASE 


ty PG ADA ste ailad i eae te w Sa estes ies w Giles Snel Sin ea ee AES Ot rear Baa as 
THREE-DIMENSIONAL INFLUENCE VALUES OF DISTURBANCE 


DISTRIBUTION OF TEMPERATURE FROM HEATED 
STRUCTURES ON PERMAFROST 


ANTM CAS ait eMUILERT SH separ tyeyset crise rtomen tenets eer telersa oh otto die leytotuerian neh <eMeLee Snore racine: 


HIGHWAY DESIGNS TO RESIST SUBGRADE MOISTURE VARIATIONS 
Be TAllan Habiburton' ach ons, ces cee si te esis igice cer auiee “ey ape: Giceieeie uielw ine) aru ators! en tere 


PREFABRICATED SUBSURFACE DRAINS 
Kent;A. Healy and-Richard.P. Long ....0008 6456 h% 6 8t ee he os ewe Be 


GRAVITY FLOW TO EXCAVATIONS AND DRAINAGE TRENCHES 
IN LAYERED AQUIFERS 


Francis G. McLean and Raymond J. Krizek........... 00. eeeeevaee 


STRESS-STRAIN BEHAVIOR OF FROZEN FINE-GRAINED SOILS 
Waddah Akili, The College of Petroleum and Minerals, Dhahran, Saudi Arabia 


The objective of this experimental study was to investigate the effect of 
below-freezing temperature, constant axial deformation rates, and soil 
type on the stress-strain behavior and strength of frozen fine-grained 
soils. Two soils were selected for this investigation: a highly plastic clay 
and a clayey silt. Samples were cored out of statically compacted soil 
cakes and were quickly frozen. Average molding densities and average 
molding water contents of test samples fell on the wet side of optimum 
conditions determined by standard Proctor curves of compaction. Constant 
axial deformation rate tests were carried out on frozen samples at -1, -5, 
-9, and -22 C, and at different constant axial deformation rates. Two types 
of stress-strain behavior were exhibited; the brittle type was associated 
with the clayey silt, and the plastic type was associated with the clay. Re- 
sults also show a strong dependency of ultimate strength (peak strength) 
derived from stress-strain curves on temperature for both soils tested. 
Ultimate strength is also shown to depend on deformation rate. The infer- 
ence may be drawn that the amount of liquid water present as a thin film 
between solid and ice surfaces and the ratio of liquid water to ice in a fro- 
zen soil are responsible for ice cementation bonds. These bonds in turn 
control the stress-strain behavior and the ultimate strength of frozen soils. 


eIT HAS BEEN SHOWN that engineering properties, such as creep behavior, strength, 
and thermal properties, of frozen soils are temperature dependent (1, 2, 3, 4, 5). When 
temperature dips below freezing, the phase composition of water in fine-grained soils 
changes accordingly. Part of the available water turns into ice, whereas the rest re- 
mains as supercooled water. 

The amount of unfrozen water in a frozen soil, its nature, and its equilibrium side 
by side with ice depend primarily on temperature, mineralogy and particle gradation, 
water content, and molding conditions of the wet soil. Different soils have different 
phase compositions at a given subzero temperature, and this composition changes ap- 
preciably with the lowering of temperature (6). The ratio of frozen to unfrozen water 
appears to depend on temperature history, salt concentration in the water (2), and 
lowest temperature reached during freezing (6). 

Many theories have been formulated regarding the mechanisms responsible for the 
presence of unfrozen water in frozen soils. Williams (7, 8) states that capillarity and 
suction properties of the soil are the cause of the unfrozen water, Others (9) explain it 
in terms of the oriented water structure. The least structured water freezes first. 
Additional cooling freezes lesser amounts of water. A film of water, although very thin, 
exists at the ice-particle interface at very low temperatures. This film, in particular, 
and the unfrozen water, in general, seem to have significant effect on stress-strain be- 
havior and strain-time behavior of frozen soils. Vialov (5) postulates that strength and 
deformation of frozen soils are controlled by cohesive bonds resulting from cementa- 
tion by ice. Such cementing is the result of the bonds between the ice crystals and the 
mineral particles that are separated by a film of unfrozen water. This type of bond is 
very unstable because it changes with any variation in the temperature field. Under the 
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influence of a load, the ice cementation TABLE 1 
bonds break up, and partial melting of ice | CHARACTERISTICS OF 2 SOILS 
occurs. The water film, which increases 


Sault 


because of phase change of ice, moves Characteristic Ste. Marie Raven i 
from zones of high stress to zones of Clay 

smaller stress where it refreezes in  jiquia limit 60 24 
time. Plastic limit 25 17 


. . n 7 Plasticity index 35 7 
To provide information on the stress Gnanitic gravity 2.79 2.69 


strain behavior of frozen fine-grained  Gradation, percentage finer 


soils required that constant axialdefor- = 7 ™ . 
mation rate tests be carried out on iden- 0.002 mm 60 8 


tically molded samples of a plastic clay 
soil and a clayey silt soil at several sub- 
zero temperatures and at different defor- 
mation rates. The test results indicate 
that stress-strain curves of the frozen samples differ in shape depending on soil type, 
test temperature, and axial deformation rate. The data show convincingly that ultimate 
strength derived from strain curves depends strongly on temperature and moderately 
on axial deformation rate. 


EXPERIMENTAL PROGRAM 
Soils Studied 


The 2 soils selected for this investigation were a red plastic clay obtained from 
Sault Ste. Marie, Michigan, and a tan glacial clayey silt obtained from a location near 
Lafayette, Indiana. The Sault Ste. Marie clay was used previously by the author (1). 
The characteristics of both soils are given in Table 1. The mineral contents of the 
clay are illite, vermiculite, and chlorite, Illite and quartz are predominant in the 
clayey silt. The compaction moisture-density relations obtained by using the standard 
Proctor test are shown in Figure 1. 

The major difference between the soils selected is the degree of plasticity and 
particle gradation. When mixed with sufficient water, the clay soil behaves like putty, 
whereas the clayey silt has no apparent plasticity when worked by hand. 


Molding of Test Samples 


The clay and the clayey silt were allowed to air-dry. Then they were processed by 
crushing and sieving until they passed the No. 40, sieve. Distilled water was added to 
the air-dried soil in predetermined 
quantities according to the desired 
molding conditions. Vigorous peri- 
By — MOED SIRT SAMPEES odic hand mixing of the moistened 
ere Eee ee soil ensured even water content dis- 
tribution prior to molding. The pre- 
pared soil was statically compacted 
in a 12-in. diameter mold specially 
par rea cancel made for this investigation. 
niet In the molding process, a prede- 
90 a termined weight of the soil-water 

| mix was placed evenly in the large 
cylindrical mold (12 in. in diameter 


0 


DRY DENSITY, Ibstt? 


80 SLAY and 6 in, high) and was subjected to 
MAX. By =89.21b/f1° static load by using a Tinus-Olsen 
LL | OM.C.= 29 % | loading machine. The desirable 
10 is 20 25 30 35 axial load was applied to the wet 
MOISTURE “CONTENT: soil in the mold through a plunger 


and was maintained constant for a 
Figure 1. Moisture-density relations determined by using the period of about 10 min until desir- 
standard Proctor compaction test. able density was reached. The 


compacted soil cake was taken out of the 
mold, sliced into 21 pieces, wrapped in- 
dividually in aluminum foil, and covered 
with wax pending later use. Cylindrical 
test samples were prepared as needed 
from the waxed pieces by punching a cy- 
lindrical mold through covered pieces 
with a hand press (Fig. 2). The cylindri- 
cal soil cores were taken out of the split 
mold (1.3 in. in diameter) and trimmed to 
the desired height. This coring process Figure 2. Compacted soil piece in wax (left), cylin- 
was selected over trimming with a soil drical split mold with soil core in it (center), and 
lathe because of its ease in producing trimmed sample in a rubber membrane with Lucite 
samples with identical cross sections, disks on-top and bottom (right). 

which are difficult to obtain when the 

trimming method is used. 

The method described in the foregoing permitted preparation of good cylindrical 
samples with no surface cracks or cavities. It is believed that the structure of the 
molded samples is the flocculated type associated with an edge-to-face particle ar- 
rangement as a result of the static method of compaction applied. Molded samples fell 
on the wet side of the standard Proctor curve of compaction (Fig. 1) for both soils. 
Dry densities and water contents of individual samples were within 2 percent of the 
average density and average water content of all samples used in this investigation. 


Sample Preparation and Freezing 


The molded soil samples were weighed and measured, Their initial water content 
was determined. Lucite disks were placed on each end of the cylindrical samples with 
specially prepared friction reducers to minimize end friction during deformation of 
samples. The friction reducers consisted of a perforated sheet of aluminum foil coated 
with a thin layer of silicone grease and covered with a polyethelene sheet on the top and 
bottom. Rubber membranes were placed over samples with tightly fitting rubber bands 
placed around the Lucite disks, thereby preventing any loss in moisture prior to or 
during testing. Samples were then placed in a pan filled with water and were stored at 
ambient temperature for a period of about 2 weeks before freezing to minimize thixo- 
tropic effects. 

Test samples were mounted in a triaxial cell, secured in place on the pedestal of 
the cell, and quickly frozen by filling the cell with a coolant. The coolant, a mixture of 
ethylene glycol and water, circulated around the cell at a constant temperature by 
pumping the coolant through a low-temperature bath. All test samples were cooled 3 C 
below test temperature, then warmed up to test temperature, and left to reach equilib- 

rium for at least 12 hours before they were 
tested. Freezing of cylindrical samples begins 


ao NGLBED DRY DENS TVS. Oo BETAS on the surface and progresses inward, Water 
z MOLDED WATER CONTENT- 33.2% in large spaces freezes first, and the ice-water 
Higa en ree Eee STUB ATION | U00 % interface will extend into smaller spaces and 
& 30 E . 
ie Ets channels until faced with too small a space to 
z \ : . 
8 vs allow ice to grow. There the water in the space 
fs aS cootine will remain supercooled as long as the temper- 
5 ‘AS ‘ ature of the sample is constant. 
i NRPS A sample warmed from a colder temperature 
2 10 =f z to a test temperature will have more moisture 
= WARMING . . 
= frozen than a duplicate sample that is cooled to 
test temperature. Figure 3 shows the variation 
oe zr ss ri ar =>. Of unfrozen water content with temperature for 
TEMPERATURE, °C a typical molded sample of Sault Ste. Marie 


clay. This was determined by the calorimetric 
Figure 3. Unfrozen water content versus method similar to the procedure used by 
temperature for Sault Ste. Marie clay. Lovell (6). Reproducibility of calorimetric 
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results for the clayey silt, used in this 
investigation, was poor and therefore is 
not reported here. 


Equipment and Testing 


The triaxial cell, with the sample 
mounted in it, was immersed in a circu- 
lar tank filled with a mixture of ethylene 
glycol and water (Fig. 4), which was main- 
tained at a constant temperature by cir- 
culating it through an adjacent cold tem- 
perature bath. Temperature fluctuations 
near sample were less than +0.05 C be- 
cause of delayed temperature response in 
the triaxial cell as compared to larger 
variations in the cold temperature bath. 
Constant axial deformation rates were 
applied to the loading ram of the triaxial 
cell by using a variable speed mechanical 
loading system. Axial loads were mea- 
sured with a load cell mounted on top of 
the loading ram and connected to an auto- 
matic recorder. Loads were recorded 
with an accuracy close to +2 lb. Uniaxial 
deformations were determined by means 
of a linear variable differential trans- 
former connected to a strip chart re- 
corder. LVDT readings were checked 
occasionally with a 0.0001 in. per division dial gage mounted on the triaxial cell and 
were substituted for chart readings whenever power fluctuations in line affected the 
accuracy of such readings. The recorder used a 10-in. wide strip of paper that per- 
mitted accurate measurement of sample deformation and served as a permanent record 
of the test. 


Figure 4. Triaxial cell in tank with coolant 
circulating. 


EXPERIMENTAL RESULTS 


Constant axial deformation rate tests were carried out on frozen cylindrical soil 
samples of clay and clayey silt at several temperatures below freezing. Deformation 
rates were maintained constant by means of a variable speed mechanical loading system. 
A summary of the characteristics of the 2 soils used is given in Table 1. Results of the 
constant axial deformation rate tests are presented in 2 groups according to the soil 


type. 
Clay Samples 


Typical stress-strain curves 
for identically prepared and 
frozen Sault Ste. Marie clay 
SAMPLE 3-2 Siete Sey samples deformed at constant 
esi ib [ &=1,55x10 in/min, deformation rates are shown in 

Figures 5,6, 7, and 8. The 
average molded dry density and 
average molding water content 
were 90.6 lb/ft? and 33.2 percent 
respectively. Test temperatures 


AXIAL STRESS, PSI 


TRUE AXIAL ain, % were -1C (Fig. 5), -5C (Fig. 6), 
-9 C (Fig. 7), and -22 C (Fig. 8). 
Figure 5. Typical stress-strain curves at —-1C for Sault Ste. At each temperature 2 tests are 


Marie clay. shown, one for deformation at a 


AXIAL STRESS, PSI 


SAMPLE 3-16 
€= 6.18 x 102in/min. 


“14 
103 in/min. 


TRUE AXIAL STRAIN, % 


Figure 6. Typical stress-strain curves at —5 C 


for Sault Ste. Marie clay. 


fast rate and one for deformation at a 
relatively slow rate, except that for tem- 
perature -9 C (Fig. 7) 3 different rates of 


deformation are shown. 


Tests shown in 


Figure 9 were conducted at approximately 
the same rate of deformation. 
ferent stress-strain curves (Fig. 9) re- 
flect the effect of temperature alone on 
the stress-strain behavior. 
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€=6.00x 10" in/min. 


SAMPLE 7-8 
E=2.02 x10 in/min. 
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Figure 8. Typical stress-strain curves at -22C for 


Sault Ste. Marie clay. 
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Figure 7. Typical stress-strain curves at -9C for 
Sault Ste. Marie clay. 
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Figure 9. Stress-strain curves at different subzero 
temperatures and approximately constant axial defor- 


mation rate for Sault Ste. Marie clay. 


Data shown in Figures 5, 6, 7, 8, and 9 reveal the following: 


1. The stress-strain behavior differs depending on the deformation rate. The initial 
part of these curves at a particular temperature does not seem to be affected by defor- 
mation rates. This is the so-called elastic part of the stress-strain diagram that 
would correspond to expulsion of air present in the sample and to elastic changes in 
the crystal lattice of the ice and mineral particles. The major portion of this strain is 
recoverable when load is lifted. 

2. The effect of deformation rate becomes well pronounced when samples begin to 
deform plastically. This is observed in all samples at all temperatures. It is apparent 
that the faster the deformation rate is, the higher the ultimate strength reached by the 
sample will be. 

3. At relatively slow rates, samples appear to gain resistance to deformation (such 
as sample 3-14 in Fig. 6 and sample 3-8 in Fig. 7) while deforming plastically. This 
is a strain-hardening phenomenon that was not observed at -22 C. Samples that de- 
formed slowly at -1, -5, and -9 C kept flowing under load, and ultimate strength, once 
reached, was almost constant the remainder of the test duration, No fracture or sud- 
den failure resulted in the preceding samples in comparison with samples deformed at 
-22 C where sudden failure took place on all samples strained at this temperature, It 
appears that samples deformed at -9 C and higher exhibit considerable viscous flow 
especially close to the melting point of ice (-1 C), whereas samples deformed at -22 C 
were less viscous and more brittle. 


Clayey Silt Samples 


Typical stress-strain curves for identically prepared and frozen clayey silt samples 
deformed at constant deformation rates are shown in Figures 10 and 11. The average 
molded dry density and average molding water content were 107 lb/ft? and 18.6 percent 
respectively. Samples 5-10, 5-5, and 5-13 shown in Figure 10 were deformed at con- 
stant temperature (-5 C) but at varying deformation rates, indicating that the faster the 
rate of deformation was, the higher the ultimate strength was. Samples 5-13 and 5-17 
reflect the effect of different temperatures (-5 and -1 C) at identical deformation rates. 
The following observations can be made for the clayey silt soil: 


1. The elastic portion of the stress-strain curves appears to be smaller than that 
of the highly plastic clay at comparable 
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T=-5°C 
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Figure 10. Stress-strain curves at —-1 C, —5 C, and dif- 


ferent axial deformation rates for Lafayette clayey Figure 11. Typical stress-strain curves at —22C for 
silt. Lafayette clayey silt. 


temperatures and deformation rates. Deformed samples exhibit well-pronounced 
stress-strain curves with well-defined peaks. 

2. Samples deformed at rates and temperatures comparable with those of the clay 
soil samples reach peak strength at smaller strains. All samples fail by sudden rup- 
ture shortly after peak strength has been reached, except at -1 C where some viscous 
flow has been exhibited. 

3. Most of the samples tested in this category demonstrated definite failure sur- 
faces with an approximate angle of 60 deg in the major principal plane. 


DISCUSSION OF RESULTS 


Based on the experimental data presented, it seems possible to distinguish between 
the behavior of those samples that deformed and failed shortly after deformation in a 
brittle fashion and those samples that deformed rather excessively with no visible 
fracture or sudden failure, which is referred to here as plastic failure. The first type 
of deformation and failure was exhibited by the clayey silt, except at a high tempera- 
ture (-1 C), and the latter type was exhibited by the plastic clay, except when the tem- 
perature was too low (-22 C) or the deformation rate was too large (2.1 x 10-!in./min). 

When a sample fails in a brittle fashion, the inference may be drawn that the ice 
cementation bonds between ice crystals and mineral particles have been overcome and 
destroyed completely. On the other hand, in a plastic type of failure these bonds seem 
to stretch permitting mineral particles to slide past one another and to rotate from 
their initial positions so that bonds are preserved and the material remains intact. 

The extent to which particles can rotate or slide past one another is dependent on the 
amount of liquid water and on the ratio of liquid water to ice present ina sample. The 
liquid water exists as a film of variable thickness between solid and ice surfaces. The 
thickness of this film depends strongly on temperature and soil type. 

Ultimate strength values derived from stress-strain curves shown in Figures 5 
through 11 as a function of test temperature at constant rates of deformation are shown 
in Figure 12. The linear variations shown are meant to be used in a relative sense. 
Both soils exhibit strong dependency of ultimate strength on temperature. For example, 
the approximate ratio of ultimate strength at -22 C to that at -5 C was 3.6 to 1 for the 
clayey silt and 3.3 to 1 for the clay soil at the same deformation rate of 2.1 x 10~° 
in./min. 

Both soils have shown moderate defor - 
mation rate effect (Fig. 12); that is, the 
ToLEEe OAV mars Ciaateza® strength increases with increasing defor - 
MOLDED WATER CONTENT~ 33.2% mation rates. For example, the ultimate 

aay strength of the clay soil increases approxi- 
E gyq? 2.10% 10 8in/min. mately 30 percent as the axial deformation 
rate increases from 2.06 x 10-*in./min to 


2400 


2000 


oe CLAYEY SILT 6.30 x 10-2in./min. This is a rate in- 


MOLDED DRY DENSITY - 107 Ib/f1° crease of the order of 33 times. 
Egyg=8.30x10 7in Amin. 


MOLDED WATER CONTENT- 18.6 % 
1200 


CONCLUSIONS 


In this investigation, constant axial de- 
formation rate tests were carried out on 
identically molded and frozen cylindrical 
soil samples in an attempt to gain infor- 
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\Lovay mation concerning the effects of soil type, 
Egy” 208#10° in/min. subzero temperature, and axial deforma- 
tion rates on the stress-strain behavior of 
ty “e ae als “20 4 frozen fine-grained soils. From the re- 


ge Ca sults of the testing, the following conclu- 


, : ; sions have been drawn: 
Figure 12. Ultimate strength derived from  stress- 


strain curves versus subzero temperature for Sault 1. Two types of stress-strain behavior 
Ste. Marie clay and Lafayette clayey silt. have emerged—a brittle type associated 
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with low plasticity soils (clayey silt) and a plastic type of failure associated with high 
plasticity soils (clay); 

2. Low plasticity soils exhibit plastic behavior when temperature is relatively high 
(-1 C), and high plasticity soils exhibit brittle behavior when temperature is relatively 
low (-22 C) or when axial deformation rates are rather fast (0.21 in./min) or when both 
occur; 

3. Ultimate strength derived from stress-strain curves depends strongly on tem- 
perature, and, based on molding conditions and soils used here, seems to increase 
about 3-1/2 times as the temperature decreases from -5 C to -22 C; 

4. Ultimate strength depends moderately on deformation rate and increases by 
30 percent when deformation rate increases from about 0.002 to 0.063 in./min; and 

5. Data in this investigation suggest that it is primarily the amount of liquid water 
and the ratio of liquid water to ice present in a frozen soil matrix that control its 
stress-strain behavior and its ultimate strength. 
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FROST BEHAVIOR OF COMPACTED SOILS 


Mian-Chang Wang, University of Rhode Island; and 
Gilbert L. Roderick, University of Wisconsin 


The frost behavior of compacted soils subjected to a unidirectional pene- 
tration of freezing temperature from the top under a temperature gradient 
of 20 F/6 in. was studied. The results indicated that rate and maximum 
depth of freezing temperature penetration, water content increase, frost 
heave, and development of ice segregation were dependent greatly on soil 
texture compositions and compaction conditions including initial degree of 
saturation, molding moisture content, and dry density. Although sand size 
content was constant, increasing the clay size content (smaller than 0.002 
mm) decreased frost heave, average water content increase, size of ice 
lenses, and depth to ice front. In general, decreasing the dry density only 
decreased the rate of frost penetration and increased the amount of water 
content increase and total heave. However, frost heave of Providence silt 
compacted by using standard AASHO compaction increased as molding 
moisture content increased regardless of the change in dry density. Frost 
heave increased linearly with average water content increase. The rate of 
increase, for a constant degree of saturation, increased as percentage of 
fines decreased. 


eFREEZING causes volume to increase in soils because of volume expansion of inter- 
stitial water and formation of ice lenses in the soils. Frost heaving in soils may result 
in engineering problems and failures in hydraulic structures, highway pavements, and 
buildings. Because of the detrimental effects of frost action, considerable progress 
has been made in clarifying the mechanism of frost action in soils (1-7), in understand - 
ing the factors affecting frost action (8-13), and in developing techniques for handling 
frost-susceptible soils (13-16). 

Frost action is a dynamic unsteady-state process, the nature of which is complex 
and involves a number of variables that are strongly interrelated. Because the in- 
tensity of frost action depends on so many factors that are indeterminate, it is almost 
impossible to predict the amount of frost heave in a soil under field conditions by 
means of laboratory tests and available theory. However, the relative potential for 
frost heave to occur in any soil could be estimated by studying the various factors that 
aid in its development. 

This study was, therefore, undertaken to investigate the effect of a number of fac- 
tors, such as molding moisture content, initial degree of saturation, dry density, and 
soil texture, on the frost behavior of compacted soils. The frost behavior that was 
studied included depth and rate of frost penetration, moisture migration due to freezing, 
frost heaving, and development of ice segregation. 


APPARATUS, MATERIALS, AND TEST PROCEDURE 


Tests were conducted by measuring temperature, overall heave, and final moisture 
content and by observing the development of ice segregation in test specimens subjected 
to a unidirectional penetration of freezing temperature under a temperature gradient 
of 25 F at the top and 45 F at the bottom of the 6-in. high test specimens. The testing 
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3.375 IN. 


3.063 IN. 
a 


apparatus included a cold chest, a gal- 
vanometer and thermocouples, dial 
gages, and an X-ray machine. A de- 
tailed description of the apparatus is 
given in the Appendix; a schematic view 
of the test setup is shown in Figure 1. 

Five soil specimens were tested si- xe oe i MICA 
multaneously. Overall heave and final Zins THERMOCOUPLES 
moisture content were determined for SPECIMEN 4 
all test specimens. Soil temperature 
readings were taken only from the spec - 
imen located at the center of the group, 
and observation of ice segregation was 
made only for the specimen located at 
the front left side. 

Three test soils were used: Provi- 
dence silt and 2 mixtures of Providence 
silt with different amounts of commer - 
cial Ca-montmorillonite plus sand-sized 
particles. One mixture had 21 percent 
Ca-montmorillonite plus 4 percent sand, 
and the other had 57 percent Ca-montmorillonite plus 10 percent sand by weight. The 
index properties of the 3 test soils are given in Table 1; the distribution of grain sizes 
are shown in Figure 2. 

Test specimens were compacted by using a drop-weight compaction machine manu- 
factured by Soiltest, nc. The device had a drop weight of 15 lb for striking the top 
piston of the compaction mold. The drop height for the weight was controlled by a clip 
arrangement on the shaft guide. 

The compaction molds, having inside diameters of 2.875 in. at the bottom and 3.063 
in. at the top and a height of 6 in., were made of 2-5/6 in. (inside diameter) by 1/4 in. 
thick Lucite cylinders. The reason for using inside-tapered cylinders was to reduce 
wall friction or resistance to heaving. Damage due to compaction was prevented by 
confining the Lucite cylinder in a steel mold having an inside diameter exactly equal 
to the outside dimension of the Lucite cylinder. 

The test specimens were compacted in 3 equal layers. Without the soil specimen 
being extruded from the Lucite mold, the specimen was capped with aluminum foil to 
prevent moisture evaporation and was then installed in the test chamber for at least 
24 hours. Temperature in the test chamber at this stage was controlled at 45 F for 
the purpose of simulating field ground temperature. 

Immediately following installation of the thermal insulation system and the heave- 
and temperature-measuring systems, testing was started by regulating the temperature 
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Figure 1. Schematic view of test setup. 


TABLE 1 
INDEX PROPERTIES OF TEST SOILS 


Texture 
: eee 
Test Specific Liquid ary Plasticity Composition 
Number Soil Gravity’ Time (4) Index Clay silt Sand 
(#)  (4#)— (#) 
1 Providence silt 2,75 28,0 24,0 4,0 9 54 37 
2 Providence silt + 21% 
Ca-montmorillonite + 4% 
sand 2.74 35.2 24.5 10.7 13 50 37 
3 Providence silt + 57% 
Ca-montmorillonite + 10% 
sand 2.71 46.8 26.0 20.8 23 40 37 


4According to International Society of Soil Science Classification System: clay size, < 0.002 mm; silt size, 0.002 to 0.02 mm; sand 
size, 0.02 to 2.0 mm. 


sib 


100 


X — PROVIDENCE SILT 


80 A— PROVIDENCE SILT + 21% Ca— 
MONTMORILLONITE + 4% SAND 


O— PROVIDENCE SILT + 57% Ca— 
MONTMORILLONITE + 10% SAND 


40 


20 


PERCENTAGE FINER THAN DIAMETER (%) 


0.5 0.2 O.1 0.05 0-02 0.01 0.005 0.002 0.001 0.0005 
GRAIN DIAMETER (MM.) 


Figure 2. Distribution of grain sizes in soils investigated. 


at the top of the test specimens to 25 F and by keeping the temperature at the bottom 
of test specimens constant at 45 F. Both 25 F and 45 F temperatures were used re- 
spectively at the top and the bottom of the test specimens throughout all tests. During 
testing, soil temperature readings were taken at total elapsed time intervals, such as 
1, 2, 3, 5, 7, 15, and 30 hours, whereas overall heave determination and ice segrega- 
tion observations were made every 24 hours. All tests were terminated after 5 days. 


EXPERIMENTAL RESULTS AND DISCUSSION 
Soil Temperature 


Because of the unidirectional penetration of subfreezing temperature from the top 
of the test specimens, soil temperature at a certain depth decreases with increasing 
time of freezing. Typical curves for the variation of the soil temperature at different 
depths with freezing time for the compacted Providence silt in both closed and open 
systems are shown in Figures 3 and 4 respectively. It is seen that the slopes of the 
curves decrease as freezing time increases, implying that the soil temperature at any 
depth drops at a rate that decreases as freezing time increases because of the de- 
creasing temperature gradient associated with increasing freezing time. In addition, 
because the curves are nearly parallel to each other, the rate of soil temperature 
change is nearly independent of the depth. The supply of excess water from the bottom 
of the test specimens in the open system causes a slower rate of cooling in the open 
system than in the closed system, as shown in Figure 5. 

Effect of dry density on the rate of freezing temperature penetration, under a degree 
of saturation of 60 percent, is shown in Figure 6. Figure 6 also shows the following: 


1. The effect of dry density was insignificant in the closed system; in the open sys- 
tem, the higher the dry density was, the deeper the freezing temperature penetration 
was. 

2. At early stages of freezing, the rate of freezing temperature penetration in the 
open system was nearly a constant regardless of the magnitude of dry density; the rate 
in the open system was slower than that in the closed system because of the effect of 
additional water supply. 


During the freezing process, the only sources of heat supply were moisture trans- 
port and the latent heat of fusion of the water in the soil. The moisture transport is 
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Figure 3. Variation of temperature at different depths with freezing time for Providence 
silt in closed system. 
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Figure 4. Variation of temperature at different depths with freezing time for Providence 
silt in open system. 
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Figure 5. Soil temperature distribution with depth. 


caused by the pressure deficiency at the ice front. The pressure deficiency is a mea- 
sure of the matrix and osmotic potentials. Because the variation of osmotic potential 
for the conditions under study is generally considered to be negligible (17), the pres- 
sure deficiency reflects essentially the matrix potential. The matrix potential is a 
function of water content, and change in water content is directly influenced by the rate 
of ice lens growth that, in turn, is controlled by the external temperature gradient and 


the physical properties of the soil. For a given external temperature gradient, 
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Figure 6. Penetration of 32F temperature as a function of freezing time for a constant degree of 
saturation in both open and closed systems. 
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therefore, the pressure gradient would be directly related to the physical properties 
of the soil including the composition of the soil and the state of the soil mass with re- 
gard to soil moisture content, density, and the like. According to Darcy's law, the 
amount of water flow is directly proportional to the soil permeability and the pressure 
gradient. It follows that the soil property is the only essential factor controlling the 
amount of moisture transfer, as long as the external temperature gradient is a con- 
stant. For a constant degree of saturation, a lower dry density results in a higher 
void ratio and leads to a higher amount of moisture transfer. Consequently, a lower 
dry density gives a larger amount of moisture transfer to raise the soil temperature 
and impede the penetration of freezing temperature. 

Latent heat depends only on the amount of water in a unit volume of soil. The lower 
the dry density is, under a constant degree of saturation, the higher the amount of water 
in a unit volume of soil will be, and, accordingly, the higher the latent heat will be. 
From the combination of these 2 effects (i.e., a large amount of moisture transport and 
a high value of the latent heat of fusion), it is therefore obvious that the depth to which 
the freezing temperature can penetrate is less for low dry density. 

The negligible effect of dry density on the rate of freezing temperature penetration 
in the closed system, as indicated by the experimental results, probably illustrates 
that the effect of the latent heat of fusion on the rate of temperature penetration is rela- 
tively less significant than the effect of excess water supply, at least for the conditions 
investigated. Temperature increase due to the dissipation of the latent heat of fusion 
of soil moisture for the range of dry densities studied was computed. It was found that 
the difference in temperature increase between the largest and smallest dry densities 
(i.e., 108 pef and 92 pcf) was not more than 1 F. Therefore, the results did not show 
appreciably the effect of dry density on the rate of freezing temperature penetration 
in the closed system. 

Figure 7 shows the effect of degree of saturation on the rate of freezing temperature 
penetration. It is seen that, under a constant dry density, the rate of freezing tempera- 
ture penetration increased as the degree of saturation decreased regardless of the 
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Figure 7. Penetration of 32 F temperature as a function of freezing time for a constant dry density in both 
open and closed systems. 
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nature of the testing system. This is as would be expected because a high degree of 
saturation, under a constant dry density, gives high moisture content that increases 

not only the amount of latent heat of fusion but also the soil permeability. Both of these 
factors favor raising soil temperature and impeding freezing temperature penetration. 
It is also seen that, with all other factors being equal, the depth to which the freezing 
temperature penetrates within a constant freezing time was greater in the closed sys- 
tem than in the open system. 

The effect of soil texture on the rate of freezing temperature penetration is shown 
in Figures 8 and 9. Figure 8 shows the variation of temperature at midheight of the 
test specimens with time. It is seen that the rate of temperature variation with time 
approached zero sooner for the soil sample containing a higher percentage of the frac- 
tion smaller than 0.002 mm. Figure 9 shows that, among the 3 soils tested, the greatest 
depth to which the freezing temperature could penetrate decreased as content of fine 
fraction increased. In addition, the rate of frost penetration increased as the percent- 
age of fine fraction increased, at least to a depth of 2.25 in. Below a depth of 2.25 in., 
however, the rate tended to decrease as the content of fines increased. 

The percentage of fine fraction in a soil mass is an important consideration because 
the pore size decreases as this percentage increases. Under a constant degree of 
saturation and dry density, the consequence of reduction in pore size is to lower the 
freezing point of the pore water and to decrease the amount of pore water that can 
move freely without the influence of particle surface forces. Asa result, the soil with 
higher content of fine fraction would be subject to a rapid frost penetration. The reason 
for the experimental results shown in Figure 9 is not fully understood, but the slower 
rate of frost penetration at a depth greater than 2.25 in. for the soil containing a higher 
percentage of fines might be due to retardation resulting from some factors such as 
excess water supply from the bottom of the test specimens. The results of the final 
water content measurements shown in Figure 10 could be used to support, at least 
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Figure 8. Variation of temperature at midheight of test specimens with time for different soil 
textures in open system. 
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Figure 9. Penetration of 32 F temperature as a function of freezing time for different soil textures 
in open system. 
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Figure 10. Final water content distribution with depth for 
different soil textures in open system. 
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partly, the reason. Figure 10 shows that the distribution of final water contents was 
such that the maximum height supplied by the excess water decreased as the content of 
fine particles increased. 


Final Water Content 


The water content distribution of test specimens, after having been frozen for 5 days, 
was determined by splitting the specimen into 6 horizontal segments of equal length. 

In general, the final water content distribution for compacted Providence silt decreased 
as depth below the surface increased, as shown in Figure 11. Because there was 
neither moisture evaporation from the top nor water supply from the bottom in the 
closed system, the average final water content should be equal to the initial water con- 
tent; that is, the moisture increase in the upper portion should be equal to the moisture 
decrease in the lower portion of the test specimen. In the open system, because of the 
excess water supply from the bottom, the average final water content was always larger 
than the initial water content. In addition, the water content increase in the top 1 in. of 
the test specimen was considerably larger than the average overall increase. 

The water content increase (i.e., the difference between the final and initial water 
contents) in the top 1 in. of the test specimens varied appreciably with varying molding 
conditions. Figures 12 and 13 show that the water content increase at the top was 
greater as molding water content increased but smaller as dry density increased re- 
gardless of the nature of the test system. The effect of degree of saturation on the 
water content increase at the top appeared to be less significant in the closed system 
than in the open system. The reason might be that, as has been pointed out before, the 
rate of moisture migration is directly proportional to the permeability of the soil, and 
the soil permeability is closely related to the degree of saturation. In the closed sys- 
tem, the moisture increase at the top resulted entirely from moisture migrating from 
the bottom portion. Because most of the moisture transportable under a given external 
temperature gradient had migrated by the time the test was completed, the effect of 
degree of saturation on the rate of moisture migration could not be reflected appreci- 
ably on the test results. In the open system, the portion of moisture increase due to 
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Figure 11. Typical final water content distributions with depth. 
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Figure 12. Influence of initial water content on water content increase in top 1 in. of test 
specimens for Providence silt. 
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excess water supply from the bottom varied with various degrees of saturation and 
resulted in different amounts of moisture increase at the end of the test. However, no 
definite trend was obtained to indicate the consequent effect of the degree of saturation. 

The general feature of final water content distribution with depth changed signifi- 
cantly with changing soil textures. It was found that the depth at which the maximum 
final water content occurred increased as the percentage of fine fraction smaller than 
clay size (0.002 mm) increased. Typical curves of the distribution are shown in 
Figure 10. The difference may be attributed to the effect of pore size. Because de- 
creasing pore size increases the energy required to cause the same amount of water 
flow, the height at which the maximum moisture content occurs increases as pore size 
increases, that is, as fine content decreases. Following the same reasoning, we can 
conclude that the average overall water content increase would be less as fine-grain 
content increased. 


Frost Heave 


Test results show that considerably less heaving developed in the closed system 
than in the open system. Some specimens in the closed system even showed shrinkage 
or contraction that resulted from consolidation within the lower portion of the test 
specimens because of soil moisture that migrated to the upper portion. 

Variation of heaving with freezing time was such that the slope of the curves de- 
creased rapidly as duration of freezing increased and eventually approached almost a 
constant value. The maximum slopes, which were equivalent to the maximum rate of 
heave, increased as dry density increased, as shown in Figure 14. 

An increase in the rate of 
frost heave with an increase in 
dry density, however, does not 

necessarily imply that the over- 
TeUetnaT ORE CHEE OO, crretn all heaving would follow the 
OPEN SYSTEM same trend. Instead, Figure 15 
¥ DEGREE OF ‘SATURATION shows that, within the range of 
A 4 conditions investigated, for a 
given degree of saturation, total 
heave decreased linearly as dry 
density increased. The slope of 
the lines increased as degree of 
saturation increased. Also, for 
a given dry density, the higher 
the degree of saturation is, the 
higher the heaving rate will be. 
Figure 15 also shows that total 
heave increased as molding 
moisture content increased for 
a constant dry density, and in- 
creased as degree of saturation 
increased. 

The effect of compaction on 
frost heaving should be a com- 
bined effect of molding moisture 
content, dry density, degree of 
saturation, and the like. A stan- 
dard AASHO compaction curve, 
shown in Figure 15, illustrates 
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Figure 15. Total heave as a function of dry density. 


that the effect of dry density was masked by the combined effect of molding water con- 
tent and degree of saturation. 

Figure 16 shows the difference in the rate of heaving for different soil textures com- 
pacted at the same dry density and degree of saturation. It is noted that increasing the 
fine grain fraction (smaller than 0.002 mm) decreased the rate of heaving. In addition, 
the intensity of heaving decreased as the percentage of fine fraction increased. 

Total heave was plotted versus average water content increase for various degrees 
of saturation in Providence silt and in various soil textures with 60 percent degree of 
saturation, as shown in Figure 17. Total heave increased almost linearly with greater 
average water content increase. The higher the degree of saturation was, for a given 
water content increase, the higher the heave was. This could be expected because the 
higher the degree of saturation was, the smaller the pore size for a constant moisture 
content was, or the larger the volume of pore water was for a constant dry density. 
The smaller pore size would give greater volume expansion under the same moisture 
increase, whereas the larger volume of pore water would make the formation of bigger 
ice lenses possible. Consequently, for a given amount of water content increase, heave 
increased as the degree of saturation increased. ; 

A marked difference in the slopes of the linear relationships for different soil tex- 
tures is noted (Fig. 17). The higher the content of fines is (smaller than 0.002 mm), 
the smaller the slope will be. The test results may emphasize the effect of the size of 
ice lenses on the overall volume change of the test specimens because, as will be seen 
later, it was observed from the X-ray radiograph that the thickness of ice layers was 
different for different contents of fine-grain fraction, which was thicker for the soil 
having less clay-size content. 


Ice Segregation 


Development of ice segregation in the test specimens was studied by using an 
X-ray radiograph and by directly observing ice lenses by splitting the specimens 
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Figure 16. Variation of total heave with freezing time for different soil textures. 


longitudinally. The X-ray radiographs were taken immediately before the test and 
every 24 hours during testing. It was found that, in general, the size of ice lenses was 
much smaller in the closed system than in the open system and that in the open system 
ice lenses decreased in size as percentage of fines increased. 

Depth to the ice front was measured directly from the radiographs. The depth in- 
creased at freezing time increased at a gradually decreasing rate, as shown in Figure 
18. Penetration of the ice front approached a maximum sooner in the closed system 
than in the open system; the maximum penetration by the end of testing was smaller in 
the closed system than in the open system. In the open system, the depth to the ice 
front increased as the degree of saturation increased while dry density was constant, 
and increased as dry density decreased while degree of saturation was constant. Depth 
to the ice front was smaller for the soil samples containing higher percentages of fines. 


SUMMARY AND CONCLUSIONS 


The frost behaviors of 3 soils compacted by using various compaction conditions 
were studied by means of a unidirectional penetration of freezing temperature under a 
temperature gradient of 20 F/6 in. (25 and 45 F at top and bottom respectively of 6-in. 
high test specimens). The results of this study have presented a more complete picture 
of the influence of factors on the frost behavior of compacted soils. It is believed that 
recognition of the importance of the many factors controlling frost action may aid in 
the selection of appropriate compaction conditions and in the improvement of methods 
in engineering practice for predicting the effects of frost action. From the results of 
this study, the following summary and conclusions appear warranted: 


Freezing Temperature Penetration 


1. The rate of freezing temperature (32 F) penetration decreased as freezing time 
increased but was faster in the closed system than in the open system. 

2. Under a constant degree of saturation, the rate and the maximum depth of freez- 
ing temperature penetration increased as dry density increased in the open system; 
however, the effect of dry density was insignificant in the closed system. 
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Figure 17. Relationship between total heave and average water content 
increase. 


3. Under a constant dry density, the rate of freezing temperature penetration in- 
creased as degree of saturation decreased. 

4, For the particular test conditions used, the rate of frost penetration increased 
above a depth of 2.25 in. as the content of grains smaller than clay size (0.002 mm) in- 
creased. Below a depth of 2.25 in., the rate tended to decrease as the content of fines 
increased. 


Final Water Content 


5. Although other factors were equal, there was a greater water content increase in 
the top 1 in. of the test specimens respectively as molding moisture content increased, 
dry density decreased, or content of fine fraction smaller than clay size decreased. 

6. The final water content distribution with depth for various soil textures was such 
that the depth at which the maximum occurred increased as percentage of fines 
increased. 


Frost Heave 


7, Increasing dry density increased the maximum rate of frost heaving but decreased 
linearly the total heave for a constant degree of saturation. The higher the degree of 
saturation was, the larger the amount of heave was. 

8. Frost heave of soils that were compacted by using a constant compaction effort 
(e.g., standard AASHO compaction) increased as molding moisture content increased no 
matter how the dry density changed. 
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Figure 18. Penetration of ice front with freezing time. 


9. Even though the content of sand size was kept constant, frost heave decreased as 
clay size content increased, i.e., as silt size content decreased. 

10. Within the conditions investigated, observations were made of linear relation- 
ships between total heave and average water content increase for Providence silt com- 
pacted at various degrees of saturation and for different soil textures under a constant 
degree of saturation. 


Ice Segregation 


11. Larger ice lenses and greater depth of ice front penetration were given by the 
open system than by the closed system. Furthermore, in the open system, the size of 
ice lenses and the thickness of ice layer decreased as percentage of fines increased. 
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APPENDIX 
DESCRIPTION OF TEST APPARATUS 
Cold Chest 


A refrigerator having an inside dimension approximately 23 in. wide by 13.5 in. long 
by 31 in. high was used to control temperature within the test chamber. The size of the 
test chamber was large enough to accommodate 5 test specimens. The test specimens, 
contained in Lucite cylinders, were separated with granular insulation, a commercial 
mica, so that freezing temperature could only penetrate unidirectionally from the top of 
the test specimens. Plates of Styrofoam were place against the door of the cold chest 
to retain the commercial mica and also to prevent both freezing temperature and out- 
side temperature from penetrating laterally into the bottom of the test chamber. 

The temperature at the bottom of the test chamber was kept uniform and constant at 
45 F by means of water circulated by a pump. The level of the circulating water was so 
controlled by an overflow cup set outside of the refrigerator that it was just in touch 
with the bottom of the test samples. Besides controlling the temperature, the circulat- 
ing water also served as a constant groundwater table as in the field condition. This 
simulated groundwater table could be used to supply additional water during freezing 
when a porous stone was set underneath the soil samples. The testing system in this 
case is called the open system. Whenever the closed system was desired (i.e., no ex- 
cess water supply), the porous stone was replaced by a thin metal plate. 


Temperature-Measuring System 


The temperature-measuring system was composed of copper-constantan thermo- 
couples and a readout unit for measuring soil temperature as well as the air and water 
temperatures in the test chamber. The readout unit was a reflecting type of galvanom- 
eter, manufactured by the Leeds and Northrup Company, No. 2436-b. The measuring 
instrument had a sensitivity of 0.1 microamperes per scale division and a system re- 
sistance of 35 ohms. The complete assembly was set beside the refrigerator so that 
readings could be taken without disturbing the test chamber temperature. 


25 


Among 5 test soil specimens, only the one at the center of the group was used to 
measure soil temperature. The thermocouples were connected to the test soil speci- 
men at intervals of 0.75 in. through the height of the specimen so that the rate of tem- 
perature penetration could be measured directly (Fig. 1). The holes on the Lucite 
cylinder for the thermocouples were staggered about the circumference of the cylinder. 


Volume Change Measuring System 


Dial gages reading to 0.001 in. were used to measure uniaxial heaving or shrinking 
of the test samples, The dial gages were fastened on horizontal rods that, in turn, were 
fixed on the wall of the test chamber. A thin metal plate was set between the dial gage 
and the test specimen to ensure a smooth surface and also to prevent the evaporation of 
soil moisture, A light, controlled by a switch located outside the refrigerator, was used 
to light the test chamber so that gage readings could be taken more easily. All readings 
could be taken through the transparent door of the refrigerator without disturbing the 
testing process. 


X-Ray Machine 


The X-ray machine was a modified medical X-ray unit manufactured by the Picker 
X-Ray Corporation. It was operated with a voltage of 80 kilovolts and a current of 
15 milliamperes, For observation of the development of ice segregation, K-ray radio- 
graphs were taken every 24 hours for the sample located at the front left side of the 
group. 


VISCOELASTIC PROPERTIES OF A LATERALLY CONFINED 
SAND-ICE SYSTEM SUBJECTED TO TEMPERATURE INCREASE 


J. T. Laba, Associate Professor of Civil Engineering, University of Windsor, 
Windsor, Ontario 


A frozen sand-ice layer, when subjected to a temperature rise, will exert 
pressure against confining boundaries. The magnitude of this pressure can 
be measured and recorded in the form of a pressure-time curve. The 
sand-ice system closely resembles a viscoelastic material, and its time- 
dependent behavior under applied pressure can be compared to the be- 
havior of the viscoelastic model. Subjecting the theoretical model to 
loading conditions equal to actual pressures exerted by a sand-ice layer 
and measured experimentally made it possible to study both the viscous 
(plastic) flow and elastic strain in a frozen soil layer. Furthermore, ap- 
plying viscoelastic model analogy to the investigated frozen sand-ice layer, 
which was radially restrained and subjected to a uniform temperature in- 
crease, made it possible to derive equations for initial modulus of elasticity 
of a frozen sand layer, effect of temperature increase on the modulus of 
elasticity, time-dependent coefficient of viscosity of a frozen sand-ice 
layer, and time-dependent functions of elastic strain and plastic (creep) 
strain. 


eTHE MAGNITUDE of the lateral thrust exerted by a frozen sand layer on rigid con- 
fining boundaries and caused by uniform temperature increase was discussed elsewhere 
by the author (1). Two different sands were used in the foregoing experimental investi- 
gation. Sand No. 1 was a crushed uniform sand from Ottawa, Illinois (uniformity coeffi- 
cient 1.5), whereas sand No, 2 was a natural well-graded sand from Paris, Ontario 
(uniformity coefficient 3.8), The grain-size distributions are shown in Figure 1. The 
soil specimens tested (4 in. high and 9-5/8 in. in diameter) were composed of sand and 
a variety of ice contents. The horizontal strain in the frozen sand was measured by 
means of BLH Type A-9 bonded resistance strain gages embedded in the sand. The test 
method used consisted of zeroing the strain bridge with the frozen soil at a constant 
initial temperature and then raising the soil temperature at a desired rate while simul- 
taneously applying and measuring the resisting radial pressure required to prevent the 
soil-ice system from expanding. 

The pressure o exerted by the frozen sand layer was found to be a function of 5 vari- 
ables: the initial temperature of the frozen layer T,, the rate of the layer's temperature 
increase 6, the time of temperature increase t, the initial porosity of the sand n, and the 
degree of ice saturation §,. 

Based on the experimentally obtained results, 3 general equations were derived for 
pressure-time curves, for values of maximum pressure developed by a sand-ice layer, 
and for period of time required by the sand-ice layer to reach its maximum pressure. 
This research study was accomplished by not introducing the viscoelastic behavior of a 
sand-ice system but by knowing only the physical soil properties n, S,, T,, and é. 

However, applying viscoelastic model analogy to the frozen sand layer and subjecting 
the theoretical model to loading conditions equal to the actual pressures exerted by the 
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sand-ice system made it possible to study TYLER! ‘STANDARD ‘SIEVE! NUMBERS 
changes occurring in the elastic and plastic 100 
strain during the time when the frozen soil 
layer is subjected to a uniform temperature 
increase and, therefore, to varying lateral 
stress. Thus, the object of the investigation 
presented in this paper is to study the time- . 70 
dependent behavior of a frozen sand layer sub- 
jected to atemperature increase but prevented 
from expanding by all-round rigid retaining 
boundaries. 
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VISCOELASTIC MODEL ANALOGY 


Typical pressure-time curves shown in 
Figures 2 and 3 clearly indicate the time- ee 
dependent behavior of the lateral pressure de- 
veloped by sand-ice specimens when they are 
subjected to temperature increase and pre- oll 
vented from expansion by retaining boundaries. 

The sand-ice resembles a viscoelastic 
material, and its behavior under applied pres- 
sure can be compared to the behavior of the 
viscoelastic model (i.e., the Maxwell model), 
which combines elastic and viscous (plastic) 
response and shows the appropriate physical 
picture of the time-dependent material. 

Let us assume that the Maxwell model is placed horizontally between 2 rigid walls 
located a distance L apart (Fig. 4) and that the model is subjected to a uniform temper- 
ature increase. Because the rigid walls cannot be moved apart (L is constant), the 
linear expansion caused by temperature increase will not take place; but, instead, 
pressure on the confining boundaries will develop. The counteractive pressure pro- 
vided by retaining walls will, in turn, apply a compressive force on the spring and the 
dashpot. Under these circumstances, the resultant horizontal linear strain must be 
zero because 
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Figure 1. Grain-size distribution. 


° 


af EXPERIMENTAL 


|| n=46.6% 
6.631 texpet/eo) | | riod 
a, exp(- ‘ 
(Specimen "z") S*50% 
——— THEORETICAL 
AL [To=O* | = 
Z J 
° 


° 


is POINTS 
= 30 s 

G $ é 

2 {| 6145 texpl-t/45) © 6= 1.413 texp(-t/25) 
+ eT oe HOP ] 6= 0.815t exp(-t/20) 
rs all ier Pail 

a / OS 

o f1|6=0.279 texpt-t/5) *Ss| | ; 

w a Sy 

c 

a 


TIME, MIN. 


Figure 2. Pressure-time curves for temperature rise of 6 F per hour and 50 
percent ice saturation from the indicated initial temperatures. 
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Figure 3. Pressure-time curves for temperature rise of 6 F per hour and 33.3 
percent ice saturation from the indicated initial temperatures. 


Ey ~ (es, Fey,) = 0 (1) 
where 
e, = strain caused by temperature increase; 
é,, = strain in elastic component, spring; and 
Sor, = strain in viscous element, dashpot. 


The analogy between the behavior of the theoretical viscoelastic model and the inves- 
tigated sand-ice specimens is obvious, if we assume that the elastic behavior of the 
frozen soil is represented by the action of the spring and the plastic flow (creep) in the 
ice matrix by the motion of the dashpot. 

Each sand-ice specimen tested (1) was circular in shape (diameter = 9-5/8 in.) and 
during the testing was subjected to a radial stress o uniformly distributed around its 
circumference. Because both the sample itself and the external loading were symmet- 
rical around the sample's center, the sample remained perfectly circular in form for 
the duration of the experiment. If we assume that the sand-ice specimen is a homoge- 
neous and isotropic material (idealization of a sand-ice system with irregular ice 
crystal orientation), then both the radial and circumferential strains must be uniform 

and equal throughout the soil sample. It fol- 
lows that each element, such as N or P shown 
in Figure 5, will be subjected to biaxial 


HEAT stresses equal in magnitude and that 
| \ 4 { j ox = oy =o (2) 
where 
ox,cy = principal stresses, and 


o = radial stress. 


According to the law of elasticity, the 
L (CONST.) value of elastic stress can be expressed as 


= (3) 
Figure 4. Viscoelastic model. feet  &§E 
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where E is the modulus of elasticity. The plastic deformation may be determined from 
Newton's law of viscosity (2) as follows: 


dew 
dt 


or 
t 
o 
en, = | 2 ue (4) 


where 7 is the coefficient of viscosity (proportionality). The strain caused by a uniform 
temperature increase is 


€y = aft (5) 
where 
a = average coefficient of thermal expansion of the material, 
6 = rate of sample's temperature increase, and 
t = time, 


Referring to the stress conditions of element N shown in Figure 5 and combining Eqs.1, 
3, 4, and 5, we have 


MG) * Fay 
aft = (l-n,) + (A1- [ dt (6) 
Ee ili jo Nt) 
where 
#, = Poisson's ratio in the elastic range; 
#2 = Poisson's ratio in the plastic range; and 
(t) = a subscript, indicating a time-dependent function. 


Because at the present time 
ly values of », and », are not avail- 
| able for frozen sand, Eq. 6 can be 
written in a different form, for the 
problem at hand, without including 
directly Poisson's ratios », and 
[ge 


STRAIN 


Ex) = modulus of elasticity, 
expressed as a ratio of 
equal biaxial stresses 
to the uniaxial strain; 
and 

= coefficient of viscosity, 
expressed as a ratio of 
equal biaxial stresses 
to the uniaxial velocity 
gradient of the plastic 

Figure 5. Compressive stresses in sand-ice specimen. flow. 
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In Eq. 7, only values of E,,, and 7,, are unknown because the values of o,, were 
obtained experimentally and recorded in the form of pressure-time curves. For each 
pressure-time curve, the relationship between pressure and time can be approximated 
by Eq. 8 (see the broken lines shown in Figs. 2 and 3). 


To, = Atexp(-t/t,) (8) 
where 
A = coefficient, psi/min; 
t = time, min; and 
t,, = time required to reach maximum pressure, min. 


From Eq. 8 it follows that 
omax = 0,3679At, (9) 
Because the values of maximum pressure o max and the required time t,, to reach 
this pressure were measured for each experiment, it was possible to calculate by using 
Eq. 10 the corresponding values of the coefficient A. 


A = o max/0.3679 t,, (10) 


Combining Eqs. 7 and 8 and differentiating with respect to t gives the following differ - 
ential equation: 


a | At exp (- t/t,,) At exp (- t/t,,) 
ad = at E + ; 
(t) Mt) 
or 
a 8 exp (t/t) E E! Bre 
a Ott, - geet et (11) 


A 


The modulus of elasticity of a sand-ice specimen, before being subjected to a tem- 
perature increase, can be obtained from Eq. 11 by introducing the following boundary 
conditions. Att = 0, = 0, and E;,, = E;; therefore, 

Ei = A/a@ (12) 


where E' is the initial modulus of elasticity of a frozen sand. 
Substituting a6 E! for A in Eq. 11 gives 
Ee) By 
Er exp(t/t,) = L=t/t,- 7 tt Tr Ft (13) 
0 (t) ct) 


Temperature Effect on E’ 


Equation 12 shows that E' is directly proportional to the coefficient A and inversely 
proportional to a and @, Because A itself isa function of T,, S,, 9, andn, it can be ex- 
pected that the modulus of elasticity E') of a sand-ice system will also be influenced by 
the previously mentioned parameters, with the possible exception of @. Values of A 
obtained from 15 experiments (sand No. 1, n = 46.6 percent, S, = 50 percent, 6 = 3, 6, 
and 9 F/hour) anda = 15 x 10-° per deg F were used to calculate coefficients of elas- 
ticity E! (Table 1). 

Figure 6 shows the relationship between E} and the initial soil temperature T, for 
3 rates of temperature increase @. It can be observed that the modulus of elasticity 
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decreases as the frozen sand's initial TABLE 1 
temperature increases. On the other ia opp sear EFFECT ON MODULUS 
hand, it appears that the rate of temper - 

ature increase 6 has very little, if any, Rate of 


: ‘ . : : Temperature Initial Coefficient Coefficient 
influence on E,. The dispersion in values Rretenca Temperature “(ci/min) 0! Elasticity 
of E! obtained for 3 rates of @, at the (deg F/min) (deg ¥) Egft IN pat) 
same initial soil temperature T,, could be ae ji ‘a6 rer 
caused by instrumental errors and also 10 0.718 9.574 
by the difference in the crystalline struc- . : — one 
ture of the ice matrix formed in the 30 0.136 1812 
otherwise identical sand-ice specimens. 0,10 0 1.631 10: B74 
The relationship between E! and T, 10 1.450 9.666 
shown in Figure 6 can be represented by = ao gen 
the following equation: 30 0.279 1.860 
0.15 0 2,718 12,080 
EL = 11.3x105- 1,030T2 in psi (14) 10 2,536 11.270 
20 1.699 7,550 
25 1.411 6,270 
30 0,544 2,418 


or less accurately ? by the straight -line Note: Sand No. 1, S; = 50 percent, n = 46.6 percent, and @ = 15 x 1075 
equation 


Ei = (11.3-0.2T,)*105 in psi (15) 


Recent literature dealing with the properties of pure ice (3, 4, 5) supports the ob- 
servations that were made on sand-ice specimens and that indicate that the elastic 
modulus of ice decreases as ice temperature increases. 


Modulus of Elasticity E(,, and Coefficient of Viscosity Net) 


An analysis of the viscoelastic behavior of the sand-ice system was carried out by 
using the specimen of sand No. 1, having n = 46.6 percent, S, = 50 percent, T, = OF, 
and being subjected to a temperature rise of @ = 0.1 F/min. The selected sample has 
been called specimen Z, and its pressure-time curves (both experimental and theo- 
retical) are shown in Figure 2. The theoretical (idealized) pressure-time curve com- 
puted from Eq. 8 by using the experimentally obtained values of A = 1.631 andt, = 
60 min gives 


Tr) = 1.631 t exp (- t /60) in psi (16) 
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Figure 6. Initial modulus of elasticity versus initial frozen sand temperature. 
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Before being subjected to a constant rate of temperature increase of @ = 0.1 F/min, 
the selected soil specimen Z was at the initial temperature of T, = 0 F and, therefore, 
had an initial modulus of elasticity E; = 11.3 x 10° psi (Fig. 6). During the progress 
of the experiment, the frozen specimen, while uniformly increasing in temperature, 
was subjected at the same time to a constantly changing stress o,,. Because both the 
stress and the soil temperature increase A T are time dependent, it can therefore be 
expected that the modulus of elasticity 1a of the specimen Z will decrease with time 
from its maximum value E(,. = Ej at t = 0 min, in the same way as the broken curve 
shown in Figure 6. Equation 14, which indicates the effect of the temperature increase 
on the initial elastic modulus of the sand-ice system, can be used as a fair approxima- 
tion of the relationship E},. versus time when the independent variable T, is replaced 
by AT = @t. This is because Eq. 14 was derived for the sand-ice system subjected to 
similar stresses and temperature changes as the specimen being discussed. There- 
fore, the elastic modulus of the specimen Z at any time after it is subjected to a tem- 
perature increase, can be expressed as 


En) = Bo- a AT? = EO - bt? in psi (16) 
where 
E; = 11.3 x 105 psi; 
AT = (@t) = increase of soil temperature, deg F; 
a = 1,030; 
t = time, min; and 


b = (262) = 10.3 


Equation 16 can be open to criticism because it involves some degree of approxima- 
tion; therefore, 2 additional, less accurate equations for Eee are included in the analy- 
sis. This will help to evaluate the effect that the changes in function E (t) will induce 
on 7) and also on the magnitude of the plastic strain «,,. 

If for purposes of comparison one assumes that the elastic modulus decreases uni- 
formly with time, then 


Ey = Ei -s,AT = Ef - st in psi (17) 
where 
E; = 11.3 x 10° psi; 
AT = (@t) = increase in soil temperature, deg F; 
s, = 20,000; and 
s = (@s,) = 2,000. 


Or, if one assumes that modulus of elasticity is a constant, then 


Be = EE = sx nF in psi (18) 
The curves showing the 3 equations for E/,, are shown in Figure 7. The circular 
points shown in Figure 7 indicate the effect of temperature increase on the soil's initial 
modulus of elasticity, when @ = 6 F/hour. The same points are shown in Figure 6 and 
were obtained from Eq. 12. Substituting in turn Eqs. 16, 17, and 18 for E (ty into differ - 
ential Eq. 13 gives the following corresponding functions for coefficient of viscosity 7() 
in Ib-min/in.?: 


. (E, ~ bt?) t 
ao 2 19 
(t) (E, ~ bt?) exp(t/t,) abt? bay = a 
E; E! - bt? . 
a cE, > we) 2 
™t) "(Bl = st) exp (t/t,) st = 
a tf ELE = 


E! Ei - st 


where 
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Nn 
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Figure 7. Modulus of elasticity versus time (or temperature) for specimum Z. 


Tr Ft (21) 
ut) exp (t/t,) + t/t, - 1 
E; = 11.3 x 105 psi, 
b = 10.3, 
s = 2,000, and 
SAND Nel t, = 60 min 
To = O°F The curves showing the values 
Seo n = 46.6% of 74) for specimen Z obtained from 
© = 6°F/hr. Eqs. 19, 20, and 21 are shown in 
$j = 50% Figure 8. 


It can be observed that the 
curves showing the 3 different func- 
tions of 7/,,are practically identical, 
except during the first 60 min when 
the curve representing Eq. 20 indi- 
cates larger values for the coeffi- 
cient of viscosity than the remaining 
two. On the other hand, there is no 
significant difference at any time 
between the corresponding values 
of ety from either Eq. 19 or Eq. 21, 
even though each equation was de- 
rived from a different function of 
Ey». The reason for this is that 
the increase in the difference be- 
tween the 2 functions of E.,,, rep- 


EQ. 20 @ 21 


(t) 
{ resented by Eq. 16 and Eq. 18 


0 20 30 (Fig. 7), changes rather slowly 
TEMPERATURE INF during the first hour after the ex- 
periment has begun. After this 


Figure 8. Coefficient of viscosity versus time (or temperature) time, the influence of Ei on Ties 


for specimen Z. appears to be negligible. 
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This leads to the conclusion that, for all practical purposes, the coefficient of 
viscosity 7(,, of the specimen Z, during the time when the confined frozen sand sample 
is subjected to a constant rate of temperature increase from 0 to 30 F, can be expressed 
satisfactorily by Eq. 21, which is the simplest of the 3 equations introduced. 


Elastic Strain «,, and Plastic Deformation (Creep) EoL 


If the functions defining the modulus of elasticity and the coefficient of viscosity for 
the investigated specimen Z are known, the horizontal elastic strain and plastic-viscous 
deformation can be determined at any time when the specimen is subjected to a temper- 
ature increase and hence to varying lateral stress. According to the viscoelastic model 
analogy, the sum of the elastic and plastic strain, at all times in the frozen specimen, 
is equal to the strain caused by the uniform rate of warming. 

The elastic strain resulting from elastic deformations of the sand-ice system can be 
calculated for the problem at hand from the following equation: 


me) 
te, = Be (22) 
7 Fey 
If we substitute for E},, the best available function E'. = E'-pt? = 11.3 x 10° - 10.3 t’, 
(t) (t) ° 
then 
At exp (- t/t,,) t exp (- t /60) 
= = in in. /in. 23 
= Ei - bt? 6.928 x 105 - 6.315 t? — i 


or, if we assume, for purposes of comparison, that E/,, is a constant E/,, = E, = 11.3 x 
105, then 


At exp (- t/t,) t -t/60 
a ; = tem t/e) in in. /in. (24) 
E 6.928 x 108 


° 


The plastic -viscous irreversible deformation in specimen Z, resulting from the ice 
flow and possibly from displacement of sand particles, can be obtained from Eq. 25, 


dC: Se - t/60 
a, 2 | ae = | ES ESE) ae (25) 
p ' ' 
lo “Mct) 0 Tt) 
or, from the more simple equation, 
fo = whitey, = 15% 10°° 4-5, in in. /in. (26) 


Figure 9 shows the relationship between plastic and elastic strain during the period 
of 300 min when the temperature in the investigated specimen Z increased from 0 to 
30 F. It can be seen that, during the first 41 min after the radially confined sand-ice 
layer is exposed to a temperature change, the plastic component of the horizontal strain 
increases rapidly but still is smaller than the elastic component. At approximately t = 
41.5 min, both strain components are equal, and, after that time, the plastic component 
representing the accumulative creep (from t = 0 min) of the sand-ice specimen is pre- 
dominant. The results shown in Figure 9 were obtained by substituting in Eq. 26a = 
15 x 10~° per deg F for the average coefficient of thermal expansion of specimen Z. 

Because «,, is directly proportional to the applied stress oc, the function of «,,re- 
sembles in shape the pressure-time curve of specimen Z. Regardless of which of the 
2, equations introduced for <«,, were used, the maximum elastic strain occurred at time 
t = ay 60 min when the specimen was subjected to the largest radial stress (o max = 
36 psi). 

It can be noticed that all time-dependent functions, namely the elastic modulus, coef- 
ficient of viscosity, plastic strain, and elastic strain, were not extended beyond the time 
t = 300 minor T = 30 F. At this time, the theoretical pressure-time curve of specimen 
7, still indicates the existence of some radial pressure, whereas the measured pressure 
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Figure 9. Elastic and plastic strain components versus time (or temperature) for 
specimen Z. 


dropped to zero. The reason for this behavior was that, when the progressively in- 
creasing average sample temperature (measured at midheight of the sample) reached 
30 F, the top part of the 4-in. thick soil layer was experiencing some thawing of the ice 
matrix. During the ice-melting process, there was no lateral thrust exerted by the 
soil-water-ice layer; therefore, the elastic and plastic strains experienced by the sand- 
ice system ceased to exist. 

Figure 9 also shows the comparison between 2 <,,-time curves obtained from Eqs. 23 
and 24 respectively. It can be observed that the difference in the magnitude of <,, given 
by the 2 curves increases progressively with time and reaches its maximum at t = 300 
min. However, at that time, not only the difference discussed but also the values of «,, 
itself obtained from either equation are negligible when compared with the magnitude 
of the creep sustained by the frozen soil specimen. If a comparison is made at t = 300 
min, then from Eqs. 23 and 26 the following values for the elastic and plastic strain are 
obtained: ¢,, = 16 win./in. and fou = 434 yin./in. On the other hand, the corresponding 
values obtained from Eqs. 24 and 26 are «,, = 3 yin./in. and «,, = 447 yin./in. 

The resulting difference between the 2 functions of <«,, has only a small influence on 
the large magnitude of <«,, and, therefore, for all practical purposes the modulus of 
elasticity E/,, for specimen Z can be assumed to be the constant value E(,, = E) without 
the corresponding values of Net and also <«,, being seriously distorted. 


CONCLUSIONS 
Based on the preceding discussion, the following conclusions can be drawn: 


1. The elastic strain-time curve of a frozen sand layer resembles in shape its 
pressure-time curve. Maximum elastic strain and maximum pressure occur at the 
same time. After reaching this time, the elastic strain decreases rapidly and becomes 
almost negligible when compared with the plastic strain. 

2. Plastic (creep) strain in a restrained frozen sand layer commences immediately 
after the frozen sand is subjected to a temperature increase. This creep strain in- 
creases continuously with time, reaching its maximum value just before the ice matrix 
thaws. 
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3. For all practical purposes without serious error, the functions of 7/,. and ep, can 
be derived by assuming that the modulus of elasticity of a frozen sand layer subjected 
to lateral stress and to a uniform temperature increase is constant and equal to the 
initial modulus of elasticity of the frozen sand layer E}. 
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THREE-DIMENSIONAL INFLUENCE VALUES OF DISTURBANCE 
DISTRIBUTION OF TEMPERATURE FROM HEATED 
STRUCTURES ON PERMAFROST 


Alfreds R. Jumikis, College of Engineering, Rutgers University, 
The State University of New Jersey 


This paper presents temperature influence values and graphs for evalu- 
ating and plotting quickly and effectively temperature fields in permafrost 
beneath heated rectangular structures. 


eTHE OUTLINE by Lachenbruch (1) of a 3-dimensional temperature distribution in 
permafrost underneath heated buildings gives a transcendental equation (Eq. 1) for 
computing the steady-state principle of thermal disturbance T (x, y, z) at any point 

(x, y, z) in the hemispace underneath an area rectangular in plan whose dimensions are 
2a X 2b or to 1 side of the vertical axis of symmetry through the center (x = 0; y = 0) of 
the rectangle. When a constant temperature T.,,,, from a heated building is maintained 
over a rectangular area on the horizontal boundary surface of the hemispace, namely 
permafrost, the equation for calculating T in a steady state of heat flow is as follows: 


T, _- 
T(x,y,2) = po { erctn| S| oe [ee 
Zz Zz 


z? + (x+a)? + (y+b)? z?> + (x- a)? + (y +b)? 
(1) 


: acten| (x + a) (y -b) F octal (x - a) (y-b) I 


zVz2 + (x + a)? + (y - b)? zVz2 + (x — a)? + (y - b)? 


coordinates of any point in the hemispace, 
long side of rectangle, and 
short side of rectangle. 


i) 
» 
tou ou 


In addition to these 4 arctan terms, there are 8 more algebraic terms involved by 
extending Boussinesq's theory of elasticity. Tables prepared by the author according 
to the theory of elasticity (2) render greater influence values than those obtained by 
Eq. 1. 

The temperature values T are to be determined either by analysis or by means of 
graphical integration, i.e., by constructing a concentric annular grid by gnomonic pro- 
jections of meridians and parallels of a hemisphere onto its north polar tangent plane 
as suggested by Lachenbruch (1). 

In 1963, Porkhayev (3) also presented an equation similar to Eq. 1 but differently 
derived for calculating temperature distribution for a 3-dimensional problem for any 
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point in the hemispace. Obviously, Eq. 1 is very cumbersome to use in calculating and 


plotting temperature fields. 

Neither in Lachenbruch's article nor in Prokhayev's paper can any data be found 
that would facilitate determining and plotting quickly and effectively the temperature 
fields in soil, namely, permafrost, from constant surface temperature such as that 


from heated structures. 


TEMPERATURE INFLUENCE VALUES 


The author prepared tabulations (Tables 1, 2, and 3) and charts (Figs. 1, 2, 3, and 4) 
of the temperature influence values for the following reasons: 


1. Heat-conduction problems in permafrost, especially permafrost thawing (degra- 
dation) from heated structures in connection with laying of foundations in the Arctic and 
Antarctic, are of interest to the disciplines of foundation engineering and thermal soil 
mechanics; 

2. The evaluation of thawing induced by various modifications of the temperature on 
the ground surface is important in problems of engineering design and logistics; and 


vu Teonst. 


INFLUENCE VALUE T/Teons} 


0.150 


= 0.00 (0.25) 5.00 


RELATIVE DEPTH Z/B 


Figure 1. Influence values of distribution of temperature from heated structure 
structure on permafrost, T/T and A/B = 1.00. 
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38. The mathematics for calculating temperature fields should be readily available 
and accessible to engineers for effective use in the design of structures in permafrost. 


The temperature distribution influence values are prepared for a square (Table 1),a 
rectangle (Table 2), anda strip (Table 3). These tables contain influence values in the 
form of dimensionless temperature ratios T/T.,,,,, as functions of relative dimensions 
a/b and relative coordinates x/a, y/b, z/a, and z/b. (Because the computer prints in 
capital letters, the tabular designations are given in capital letters, such as A/B, X/A, 
Y/B, Z/A, and Z/B.) Data given in the tables pertain to temperature fields in the ver- 
tical x = 0 planes that pass through the center (x = 0, y = 0) of the rectangles. The 
tables are prepared for a constant surface temperature of T.,,,,=1C. For tempera- 
tures T.,,,,, other than 1 C, the tabular influence values must be multiplied by the tem- 
perature T in hand. 
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Figure 2. Influence values of distribution of temperature from heated 
structure on permafrost, T/T and A/B = 2.00. 
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TABLE 1 
INFLUENCE VALUES OF DISTRIBUTION OF TEMPERATURE FROM HEATED STRUCTURE ON PERMAFROST, 
= (Y/B) AND A/B = 1.00 


T/T, 
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z/8 
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0.04462 
0204299 
0004142 
0203990 


0.03844 
0003704 
0.03569 
9203441 
0003317 


0.03199 
0.03086 
0.02978 
0.02874 
0.02775 


12926°% 
*.G259r 
1092593 
202420 
0.02341 


9.81932 
f.01e73 
AeSLALT 
Cel1763 
Me M71 


Ne M1 A62 
OeM1614 
9.01568 
N.11524 
Me11482 


O,91441 
9.01402 
0.91365 
9.91328 
9291294 


201260 
9.91228 
0.01197 
02e01167 
0.91138 


9001110 
0091083 
0.91956 
9091031 
9291007 


0.90984 
0.90961 
0.90939 
0.90918 
0200897 


0.00000 


0.90742 
. 01464 
0092144 
02769 
0.03326 


0203809 
0004215 
0004547 
0.04807 
C.0s50C2 


0,05139 
0.05225 
De 95267 
0.05271 
0.05245 


Ce05193 
9205121 
0.05033 
0.04931 
0.04821 


204703 
0.04580 
0.04454 
0004327 
Oc 04200 


0.04073 
0003947 
0.03823 
0.03702 
0.03584 


0.03468 
0203356 
0.03247 
0603142 
0.03040 


202942 
0202647 
0.02756 
0.02668 
0.02583 


AOPSCL 
6902423 
M2347 
Ca 02274 
1.97205 


0.92138 
(22073 
e201 
NeOISEL 
f..11998 


Ge VWA29 
Fe1756 
01735 
PeX16 C6 
Ue M1639 


291593 
7.01549 
MeCISN7 
9.01467 
001427 


9201399 
0091353 
0.01318 
9.01285 
9.01252 


9.01221 
0.01190 
OoO11l61 
9091133 
0.91105 


0201979 
DeD1INES 
0.01029 
0.01005 
9200962 


0,00959 
0.00938 
920C917 
9.90897 
0.00877 


9.00C00 


0.00513 
2.01016 
C2 01497 
9.01948 
Ce 02363 


C.02736 
On O3C64 
02 03346 
0.03582 
0203776 


0.93928 
0.04043 
0. 04124 
0.04175 
. 04200 


0.04202 
0.04185 
0204152 
0204105 
Oc 04046 


0.03979 
0093904 
0.03824 
0.03740 
De 03653 


0.03563 
0.03473 
0.03362 
0.03291 
0.03201 


02003112 
0. 03024 
0.02938 
0.02854 
0.02772 


9202691 
Oo 02613 
Oo 02537 
9.02463 
0.02392 


N,02322 
QO. 32255 
©.02199 
ee I2127 
Me V20F7 


,020CR 
C. NOD 
9.91897 
0.91844 
Co91793 


3091743 
0.01696 
AM. 165% 
r.16n6 
Oe 71563 


9.91522 
0.01462 
0.01443 
De 1406 
01370 


9.01336 
C.013C2 
0.01270 
9.01238 
M. 91208 


0.01179 
Oo91151 
0.91123 
9.01097 
0.91071 


Co 91046 
9.91022 
0.00999 
290977 
0200958 


0.09934 
0. 00913 
0.00893 
0.00874 
9200855 


TABLE 2 
INFLUENCE VALUES OF DISTRIBUTION OF TEMPERATURE FROM HEATED STRUCTURE ON PERMAFROST, 
= f(Y/B) AND A/B = 2.00 


2(A) = FOOTING LENGTH (X—DIRECTION) 


T/T 


z/8 


Oo2 


025 


120 


lel 


900 


1000 


10.2 
1004 
1006 
1068 
11.0 


Alo2 
Lle4 
1106 
1108 
1220 


const 


1.00000 


0.92905 
0.85941 
0079226 
0.72852 
0.66883 


0.61353 
0.56273 
0051634 
0047416 
0243591 


0040126 
0.36995 
Oo 34160 
0.31595 
0.29271 


0o27164 
Oo25251 
Oo235i1 
0.21928 
0.20483 


0019164 
0.17958 
0.16852 
0.15838 
0.14906 


0.14048 
0013257 
0.12826 
Oo 11851 
0.11225 


0.10645 
0.10106 
0.09606 
0.09139 
0.08704 


0.08298 
0.07919 
0.07564 
0007231 
0206919 


0.06350 
0.95047 
0.05299 
0.04999 
0.04641 


0.04319 
0.04029 
0.03766 
0.93526 
0.03311 


0.03113 
0.02932 
0.02767 
0.02614 
0.02474 


0.02345 
0.02225 
0.02114 
0.002011 
0.01915 


0201626 
0.01743 
0.01666 
0.01593 
0.01525 


0.01461 
0001402 
0.01345 
0.01292 
0.01242 


0.01195 
0.01151 
0.01109 
0.01069 
0.01031 


0.00995 
0200961 
0.00929 
0.00898 
0.00869 


1.200000 


0. 92496 
Oc 85175 
0. 78190 
0071646 
0065601 


0.60069 
0.55039 
0.50482 
00 46363 
0042643 


0.39283 
Oe 36247 
0.33501 
0631015 
Oo 28761 


0.26716 
0024857 
0.23165 
0.021623 
0.20214 


0.18927 
0017747 
0016666 
Oo 15672 
0014758 


0013915 
0.13138 
0.12420 
0011755 
Oo11139 


Oo 10567 
0.10036 
0.09542 
0.09081 
0. 08652 


0. 08250 
0.07875 
0.07524 
0. 07194 
0.06885 


0.06322 
0.05822 
0.05378 
0.04981 
0004626 


0. 04306 
0.04017 
0003756 
0.03519 
0.03303 


0003106 
0002926 
0202761 
0202609 
0.02470 


0.02341 
0.02221 
QoO2111 
0.02008 
0.01913 


0.01624 
0.01741 
0.01664 
0.01591 
0.01524 


0.01460 
0.01400 
0201344 
0201291 
0.01241 


0001194 
0201150 
0.01108 
0.01068 
0.01030 


0.00995 
0.00961 
0.00928 
0.00898 
0.00869 


A/B= 


1200000 


0.90873 
0.82226 
Oo 74368 
0.67398 
0061274 


0055893 
Oo51147 
0. 46937 
Oo 43162 
0039816 


Oo 36787 
0234051 
0031574 
Oo 29325 
0.27280 


0025416 
0023715 
Oo 22161 
0.20738 
0019433 


0.18235 
0.17134 
Oo 16121 
Oo 15187 
0.014325 


0.13528 
0012791 
0.12108 
0011475 
0. 10886 


0.210339 
0.09830 
0.09354 
0.08911 
0.08497 


0.08109 
0.07746 
0.07405 
0. 07086 
0.06786 


0.06238 
0205751 
0.05317 
0204929 
0204580 


0204266 
02039863 
0003726 
Oo 03492 
0.03280 


0.03086 
0.02908 
0.02745 
0202595 
0. 02456 


0.02329 
0202211 
0.02101 
0.01999 
0201905 


0.01817 
0201735 
0.01658 
0.01586 
0.01518 


0001455 
0.01396 
0. 01340 
0201287 
0.01238 


0.01191 
0.01147 
0.01105 
0.01065 
0201028 


0.00992 
0.00959 
0200927 
0. 00896 
0.00867 


200 


1.00000 


0.85375 
0073515 
0064649 
0.57910 
0.52550 


0.48101 
0.44286 
0.40938 
0037954 
0.35265 


0.32625 
O02 30600 
0.28563 
0. 26695 
0.24978 


0.23397 
0.21941 
0.20598 
0.19358 
0.18213 


Ool7154 
0.16173 
0.15265 
0014423 
0.13642 


0.12917 
0.12242 
0.11614 
0211029 
0.10484 


0.09975 
0209499 
0.09054 
0.06638 
0.08247 


0. 07881 
0.07538 
0.07215 
0.06911 
0.06625 


0.06101 
0205634 
0.05217 
0.04843 
0.04506 


0.04201 
0203926 
0.03676 
Co03449 
0203241 


0203051 
0.02877 
0.02718 
0202570 
0202435 


0.02309 
0.02193 
Qo 02085 
0.01985 
0.01892 


0.01805 
0.01724 
0001648 
0.01577 
0.01510 


0001448 
0201389 
0201334 
0001281 
0.01232 


0.01186 
0001142 
0.01101 
0.01061 
0.01024 


0.00989 
0.00955 
0.00923 
0.00893 
0.00864 


0.e50000 


0047752 
0.45517 
0243310 
0041143 
0239028 


0.36974 
0. 34990 
0. 33083 
0.31258 
0.29517 


0.27863 
Oo 26296 
0.24815 
0.23420 
0022107 


0.20873 
Oel9716 
0. 18632 
O.17617 
0216667 


0.15778 
0014946 
0014169 
0213441 
0.12761 


0.12124 
0.11528 
0.10970 
0010446 
0.09956 


0.09495 
0.09063 
0.08657 
0.08276 
0.07917 


0.07579 
0.07260 
0006960 
0.06677 
0.06409 


0.05918 
0.05478 
0.05062 
0.04726 
0.04405 


0.04113 
0.03849 
0.03609 
0.03389 
0.03168 


0.03005 
0002836 
0.92680 
0.02537 
0002405 


0.02282 
0.02169 
0.02063 
0.01965 
0.01874 


0.91768 
0.01708 
0001634 
0.01564 
0.01498 


0.01437 
0.01379 
0.01324 
0601273 
0001225 


0.01179 
0601135 
0601094 
0.01056 
0.01019 


0.00984 
0.00951 
0.90919 
0.00889 
0.00860 


2¢8) 


0.00000 


0.10083 
0.17431 
0.21843 
0.24212 
0.25309 


0.25625 
0.25452 
0.24970 
0.24293 
0223495 


0.22625 
0.21719 
0.20800 
0. 19885 
Oo 18986 


0.18110 
0017264 
02 16450 
0.15672 
0014929 


0.14222 
0.13551 
0212915 
0.12313 
Oo11744 


0.11205 
Oc 10696 
0.10216 
0209762 
0209333 


0.089268 
Oo 08546 
0.08185 
0.07843 
0.007521 


0007215 
0.06926 
0.06653 
0.06394 
0.06148 


0005694 
0.05286 
0.04917 
0.04583 
0204280 


0.04004 
0003754 
0203524 
0.03315 
0.03123 


0.02946 
0.02784 
0.02634 
0.02495 
0.002367 


0.02248 
0202138 
0202035 
0.01940 
0.01851 


0001767 
0.01689 
0001616 
0.01548 
0.01484 


0.01423 
0001367 
0.01313 
0201262 
0.01215 


0.01170 
0.01127 
0.01087 
0.01048 
0201012 


0.00977 
0.00945 
0200913 
0.00884 
0.00856 


= FOOTING WIDTH (Y-DIRECTION) 


X7A= 0200 


0.00000 


9.04442 
0.08439 
Oo11714 
0.14198 
0215962 


0.17128 
0.17624 
0.168160 
0.18226 
0.18093 


0.17813 
0017429 
0.16970 
Oo 16462 
9015922 


0.15365 
0.14800 
0014236 
0.013679 
0213133 


0.12602 
0.12087 
0.11590 
Ooll112 
0.10654 


0.10215 
0.09795 
0.09395 
0209013 
0.08649 


0.08302 
0.07973 
0.07659 
0.07360 
0.07076 


0.06806 
0206549 
0.06305 
0.06072 
0.05851 


0205439 
0.05065 
0.04726 
0004417 
0.04135 


0.03877 
0.03642 
0.03426 
0.93226 
0203045 


0.02877 
0202722 
0.02578 
0202445 
0.02322 


0.02208 
0202101 
0.02002 
0.01910 
0.01623 


0.01742 
0.01667 
0201596 
0.01529 
0.01466 


0.01407 
0001352 
0.01299 
0.01250 
0.01203 


0.01159 
OcoOl117 
9201077 
0.01039 
0.01004 


0.00970 
0.00937 
0.00907 
0.00878 
9.090850 


1.75 


0.00000 


0.02541 
0.04946 
0.07104 
0.08950 
Co10458 


0.11638 
0.12517 
9013135 
0.13532 
0013747 


0.13814 
0.13764 
0.13621 
0.13407 
0.613139 


0.12830 
0.12493 
0.12136 
Ocll767 
0011391 


0.11014 
0210639 
0010269 
0.99905 
0.09550 


0209205 
0.08871 
0.08547 
0.08235 
0207934 


0.07645 
0.07368 
0.07102 
0.06846 
0.06602 


0.06367 
0.06143 
0.05929 
0.05724 
0.05527 


0.05160 
0.046823 
0.04515 
0.04233 
0203973 


0.03735 
0.03517 
0.03315 
0.03129 
0.02957 


0.02799 
0.02652 
0.002515 
0.02389 
0.02271 


0.02162 
0.02059 
0.01964 
0.01675 
0.01792 


OeO17135 
0.01640 
0.01571 
0201507 
0.01446 


0.01388 
0001334 
0.01283 
0601235 
0201189 


0.01146 
0601105 
0.01066 
0201029 
0200994 


0.00961 
0.00929 
0.00899 
0200870 
0.00843 


22.00 


0.00000 


9.01636 
0.03216 
9204690 
0.06021 
QcO7184 


0298169 
0208977 
0009618 
0.10105 
0.10456 


0.10687 
0.10817 
0210860 
OQ. 10833 
0.10746 


0.10613 
0019442 
0.10242 
0.210020 
0.09781 


0.09530 
0.09271 
0.09009 
0.08744 
O-. 08480 


0.08216 
0.0796C 
0.07707 
0.07459 
0.07217 


0. 06982 
0.96754 
0.06533 
0206320 
0206113 


0.05914 
0.005722 
0.05537 
0.05359 
0295188 


0.04865 
0294566 
0.04290 
0.94035 
0203799 


0.03582 
02 03380 
0.03194 
0.03021 
0.02861 


Oe02713 
0.02574 
0202446 
0.02326 
0.02214 


0.02110 
0.92013 
0.091922 
0.01836 
0291756 


0.91681 
0001611 
0.01544 
0201482 
0201423 


0.01367 
0.01315 
0091265 
0201216 
OoO1174 


0001131 
0.01091 
0201054 
0.91017 
0.00983 


0.00951 
0.00920 
0.0¢890 
0.00862 
0.00835 


0.00000 


92091129 
0.02231 
C293280 
0204255 
0.05139 


0.05924 
0206603 
O.O7177 
0.07648 
0.08025 


0.08313 
0.08523 
0.08662 
0208739 
0. 08764 


0.08744 
0098686 
0.08596 
0.08480 
0. 08344 


0.08190 
0.08024 
0.07848 
0.07664 
0.07476 


0.07285 
0.07092 
0. 06900 
0.06709 
0.06520 


0.06333 
0.06150 
0.05971 
0.05796 
0.05625 


0.05459 
0.05297 
0205140 
0.04989 
0004841 


0.04561 
C243C0 
0.94056 


0.03828. 


0093616 


0.03419 
0.03236 
0.03065 
0.02906 
0.02758 


0202620 
0.02491 
0.02371 
0.02258 
0.02153 


0.02054 
0.01962 
0.01875 
0001794 
0.01718 


0.01646 
0.01578 
0.01514 
0.01454 
0001397 


0.01344 
0.01293 
0.01245 
0201199 
0.01156 


0.91115 
0001077 
0.01040 
C.01005 
0.00971 


0.00939 
0.00909 
9.00880 
0.90853 
C.90826 


0200900 


9.00817 
0001619 
0.02391 
0.03122 
0293801 


0.04420 
0.04975 
0.05462 
0.5883 
0.06237 


0.06529 
0.06762 
0.06942 
0.07072 
0.207158 


0.07275 
0.07218 
0.07201 
0.07159 
0207095 


0007013 
0.06916 
0. 06806 
Oc 06686 
0.06558 


0.06425 
Oo 06287 
0.06145 
Oe 06002 
0.05858 


0.05714 
9205570 
0.05428 
0.05287 
0.05149 


0.05012 
0.04879 
0004748 
0204620 
0.04496 


0.04257 
9.94030 
0203817 
0.03617 
9.03428 


0.03252 
0.03086 
9.02931 
0.02786 
0.02650 


0.02522 
0. 92403 
0.02291 
0.02186 
0.02087 


0.01995 
0.91908 
0.01826 
0.91749 
0001676 


0.91608 
0.01543 
0. 01482 
0.01424 
0.01370 


0.01318 
0.01269 
0.01223 
0.01179 
9.011386 


0.01098 
0.01060 
0001024 
0.90990 
0.00958 


0.00927 
0.00897 
0.00869 
0.90842 
0.00817 


TABLE 3 
INFLUENCE VALUES OF DISTRIBUTION OF TEMPERATURE FROM HEATED STRUCTURE ON PERMAFROST, 


Tseng 


zZ/3 


1.00000 


293655 
0.87433 
0.81445 
0.75776 
0.70483 


0.65596 
0.61120 
0057045 
0.53348 
9.250000 


0.46971 
0044228 
0.41743 
De 39486 
0. 37433 


0.35562 
0.33851 
CQ. 32283 
0. 30843 
0.29517 


Oo 28293 
0.27160 
0.26110 
Oo 25133 
0.24224 


0.23375 
0.22581 
0.21638 
0221140 
0.20483 


0.19865 
Oo 19262 
Oe 18732 
Oo 16211 
Oo lT7IT 


0017249 
Oo 16804 
Oo 16382 
0215979 
0015596 


0.14680 
0.14227 
0013627 
13076 
0212566 


0.12095 
0.11657 
0211250 
0010869 
0010514 


0210180 
0.09867 
0.09573 
0209295 
0. 09033 


0.98786 
0.08551 
0.08329 
0.068117 
0.07916 


0.07725 
0.07543 
0.07369 
0207203 
0.07044 


0.06892 
0.06747 
0.96607 
0.06473 
0.06345 


0.06221 
0.06102 
0.05968 
9.05878 
0.05771 


0.05669 
0.05570 
0.05474 
0205382 
0.05293 


2(A) = 


1.00000 


0.93240 
De 86655 
0.80390 
0074546 
0.69171 


0.64277 
0259648 
0.55852 
0.52251 
0.49006 


0046077 
0.43430 
0.41032 
0238855 
0.36873 


0.35063 
0033407 
0.318848 
0.30490 
0.29201 


0.28009 
0.26905 
0.25880 
0024926 
0.24036 


0.23204 
0222426 
0.21696 
0.21010 
0.20364 


0.19756 
0.19162 
0.18639 
0.16125 
0.17638 


Ool 7176 
O. 16736 
0016318 
0.15920 
0.15541 


0.14833 
0.14185 
0.13590 
0.13043 
0.12537 


0.12069 
0.11633 
0.11228 
0.10850 
0.10496 


Oe 10164 
0.09853 
0.09560 
0.09283 
0.09022 


0.08775 
0.08541 
0.06320 
0.08109 
0007909 


0.07718 
0.07537 
0007363 
0.07197 
0.07039 


0006887 
0.06742 
0.06603 
0.06469 
Oo 06341 


0.06218 
0.06099 
0.05985 
0.05874 
0.05768 


0.05666 
0.05567 
0.05472 
0.05379 
0.05290 


2 f(Y/B) AND A/B = 1,000 = oo 


FOOTING LENGTH 


A/B=100020 


1.00000. 


0.91598 
O. 83669 
0.76515 
0.70227 
0.647568 


Q. 60002 
0.55845 
0.52185 
0.48939 
0046042 


0.43439 
0.41089 
0. 38957 
0.37016 
0.35242 


0033615 
0032118 
Oe 30739 
0.29463 
0.28281 


0.27183 
0.26162 
0.25209 
0024319 
0. 23485 


0.22704 
0221970 
0.21280 
0020629 
Oo 20015 


0.19435 
0.16886 
Oo 18366 
0017873 
00 17405 


Oe 16959 
0.16535 
Oe 16132 
0015746 
0.15378 


0.14690 
0.14060 
0.13480 
0.12945 
0.12450 


Oe11991 
0.11563 
Oe11165 
0.10793 
0.10444 


9210117 
0.09610 
0.09520 
0.09247 
0.08989 


0.08745 
0.08513 
0. 08294 
0.08085 
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Figure 3. Influence values of distribution of temperature from heated 
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HIGHWAY DESIGNS TO RESIST SUBGRADE MOISTURE VARIATIONS 


T. Allan Haliburton, School of Civil Engineering, Oklahoma State University 


Findings, conclusions, and recommendations of a 6-year study to measure 
Oklahoma subgrade moisture conditions are presented. General trends in 
Oklahoma subgrade moisture are described, as are the effects of highway 
components on subgrade moisture behavior. Resistance of various highway 
components to subgrade moisture effects is also discussed, and specific 
recommendations are made for highway design and construction on expan- 
sive subgrades. Two modes of behavior, subgrade moisture accumulation 
and subgrade moisture variation, were found to exist under Oklahoma pave- 
ment systems. Combinations of these 2 modes also occurred. Rigid pave- 
ments were found to be extremely susceptible to cracking fromvertical and 
lateral subgrade expansion, with resulting infiltration of water to the sub- 
grade. However, rigid pavements performed well until upper subgrade 
moisture contents approached the liquid limit. Flexible pavements were 
more resistant to cracking from vertical and lateral subgrade expansion 
but were susceptible to distress when subgrade moisture contents approached 
the equilibrium value of 1.1 to 1.3 times the plastic limit. Recommenda- 
tions include the use of a flexible and impervious base, subbase, or mem- 
brane component below the wearing surface of any pavement on expansive 
subgrades. Improved shoulders at least 8 ft in width are recommended, 
as is the establishment of equal drainage behavior on both sides of the 
highway section. More detailed subgrade soil testing techniques are also 
recommended. 


eTHIS PAPER summarizes findings and recommendations from a 6-year study of sub- 
grade moisture conditions under existing Oklahoma highways. Conducted by the School 
of Civil Engineering at Oklahoma State University, the project was formally initiated in 
June 1964, After preliminary planning (1), the first of 52 field test sites was installed 
under existing pavement in June 1966 (2). Subgrade moisture and density data were 
collected periodically with nuclear depth moisture and density probes, which were cali- 
brated for particular Oklahoma soil conditions (3, 4). In addition to subgrade moisture 
and density data, information concerning the following factors was compiled: precipita- 
tion, air temperature, soil type, highway design and construction history, pavement and 
shoulder performance, traffic, subgrade temperature, and pavement heave or settle- 
ment or both. 

Data were reduced immediately on collection, and continual evaluations were carried 
out during the 4-year data collection phase of research. Periodically, overall evalua- 


continued in June 1970. 


SUBGRADE MOISTURE BEHAVIOR IN OKLAHOMA 


This section summarizes research findings concerning subgrade moisture behavior 
observed under existing Oklahoma highways. More detailed analyses and data from 
which conclusions were drawn are available elsewhere (5, 6, 7, 8, 10, 11, 12, 13). 


Sponsored by Committee on Environmental Factors Except Frost and presented at the 50th Annual Meeting. 


45 


46 


General Conditions Existing in Oklahoma 


All sites were located on existing highways in the central and north central- 
northeastern part of Oklahoma. Annual rainfall in this region varies between 10 and 
40 in., increasing from west to east. Monthly rainfall amounts are highly seasonal; 
more than 12 in. may fall during the spring in the eastern part of the state, and only 
2 in. or less may fall during the winter at some western locations. The average monthly 
mean of air temperatures range from over 80 F during July and August to approximately 
30 F during January and February. In many regions of the state, the water table is 
located close to the surface and exhibits seasonal movement, rising during winter 
months and falling during summer months. A period of drought in Oklahoma ended in 
1965, just before installation of the first field research sites. 

Clays and clay shales often utilized as subgrade material in this portion of the state 
may vary somewhat, but their origins and stress histories are similar. These soils 
are normally classified by the AASHO system as A-6 to A-7 and by the Unified system 
in the upper portion of the CL range or in the CH range. The majority of cohesive soils 
are preconsolidated (usually by desiccation) and, after remolding, exhibit higher volume 
change potential than might be predicted from plasticity alone, particularly if compac- 
tion is dry of optimum. 


Observed Trends in Subgrade Moisture Behavior 


Data obtained from the study indicate that 2 basic types of moisture behavior, sub- 
grade moisture accumulation and subgrade moisture variation, exist in expansive 
Oklahoma subgrades. Factors responsible for each type of behavior are summarized 
in the following sections. 

Subgrade Moisture Accumulation—Subgrade moisture contents under new and rela- 
tively new construction and older existing pavements with wide improved shoulders and 
excellent pavement ratings tended, after a short initial period at construction moisture 
contents (usually below optimum compaction moisture), to increase without significant 
variation during an 18- to 24-month period until an equilibrium moisture content of 
approximately 1.1 to 1.3 times the subgrade plastic limit was reached. The word equi- 
librium is actually a misnomer because, in most cases, variations in moisture content 
begin to occur after the equilibrium point has been reached. The prime criterion for 
subgrade moisture accumulation without significant variation was found to be an im- 
pervious pavement system. 

For these impervious pavements, moisture accumulation is thought to occur prima- 
rily from capillary sources, i.e., prevention of evaporation, but some sites obtain 
moisture by infiltration from outside paved shoulders. Subgrade temperature gradients 
were found to cause only small moisture changes (5). 

Moisture accumulation noted under older construction in excellent condition probably 
resulted from the end of a drought cycle in 1965 just before the first research sites 
were installed. Increased availability of moisture from rainfall and rising water tables 
probably caused this behavior. 

At sites with less than excellent pavement ratings (usually indicative of pervious 
pavement or open joints), open shoulders, or fair-to-poor drainage, subgrade moisture 
variations were found to occur superimposed on the overall accumulation trend. The 
majority of these moisture variations were caused by rainfall infiltration and evapora- 
tion, usually from outside the shoulders or through the pavement surface, but they did 
not appear to halt the rate of moisture accumulation, which continued in most cases 
until subgrade moisture contents were above the plastic limit. After reaching this equi- 
librium value, moisture contents at the first 2 types of sites were subject to large vari- 
ations. Initial accumulation behavior noted at these sites is also thought to result from 
general drying of subgrade caused by the previously mentioned drought. 

Subgrade Moisture Variation—At most research sites where moisture accumulation 
was not in progress, moisture variations were found either to occur seasonally in annual 
cycles (with maximum moisture contents occurring during winter months) or to be pre- 
cipitation dependent. Most research sites where purely seasonal moisture variations 
occurred were on pavement that was rated as being excellent and impervious, and 
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moisture variations had little relation to measured precipitation, particularly under 
pavement centerlines, It was thought initially that these variations were temperature 
induced, but, although temperature-induced moisture migration does occur in Oklahoma 
subgrades, it is of relatively low magnitude, causing only a 1 or 2 percent engineering 
moisture content variation annually (5). In almost all cases for sites where moisture 
variations were seasonal and could not be related to rainfall, the great majority of 
moisture variation was found to be caused by seasonal water table movement, moving 
the zone of capillary rise, or else by delayed infiltration from areas adjacent to the 
pavement, caused by particular highway drainage conditions. Most seasonal variations 
under impervious systems did not exceed 5 percent engineering moisture content, and 
many were half this value. 

Seasonal trends were also noted to occur at sites located on pervious pavements, but 
cyclic variations were affected considerably by precipitation. Precipitation-dependent 
variations were also noted to occur at most sites having open shoulders, despite pave- 
ment condition. Many sites where precipitation and evaporation affected subgrade soil 
moisture were located on rigid pavement sections modified by asphaltic concrete 
overlay. 

As a general rule, upper subgrade moisture variations generally lagged rainfall by 
6 to 8 weeks, with longer times being required for variations to occur at deeper depths. 
Magnitude of variations was highly dependent on overall pavement condition, whether or 
not sealed shoulders were present, and on type of base and subbase material used in the 
pavement section. The magnitude of variations caused by precipitation and evaporation 
was almost always higher than that of variations produced by seasonal trends. At some 
sites the variations exceeded 10 to 15 percent engineering moisture content during 6- 
month periods and produced easily noticeable changes in pavement and shoulder condi- 
tion, At some sites, rainfall infiltration produced moisture contents close to the sub- 
grade liquid limit, with resulting loss of subgrade support and rapid pavement 
deterioration. 


Effect of Subgrade Moisture Behavior on Soil Volume Change 


Oklahoma cohesive soils are particularly subject to volume change as moisture con- 
tent varies. Volume change data were obtained from moisture measurements, subsur- 
face bench marks installed at research sites (7, 11), and general observations of pave- 
ment behavior. At most research sites only the upper 5 to 7 ft of the subgrade made 
any large contribution to subgrade volume change; moisture contents below this level 
remained relatively constant. Vertical movements were not extremely large. An em- 
pirical relationship of 1 in. of pavement heave for a 10 to 12 percent increase in engin- 
eering moisture content was developed from obtained data. This correlation was ex- 
trapolated from smaller recorded measurements (1/2 to 5/8 in.) obtained at several 
sites on A-6 and A-7 subgrades. Moisture contents were in the vicinity of the subgrade 
plastic limit, 

However, lateral subgrade expansion probably affects pavement performance to a 
larger extent than vertical swelling. The unit swelling potential of cohesive Oklahoma 
subgrades, compacted under normal conditions, is lower in the lateral direction than 
in the vertical (14). However, lateral subgrade expansion takes place over a 24- to 
40-ft width, and resulting movements were found, in some cases, to exceed 3 to 5 in. 
Tensile stresses, produced by lateral subgrade expansion, in the pavement system 
caused longitudinal cracking of the subbase, base, and pavement structure. This crack- 
ing was aggravated by flexural stresses from differential vertical heaving, even though 
the heaving was of rather low magnitude. 

Different vertical movement conditions were encountered under moderate and high 
fill sections. For these sections, moisture movement was downward and toward the 
outer edges of the fill, was caused by shoulder drying, and resulted in shoulder settle- 
ments of up to several inches in magnitude during dry periods. Some rebound was ob- 
served in succeeding wet periods. 
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INFLUENCE OF HIGHWAY DESIGN FACTORS 
ON SUBGRADE MOISTURE BEHAVIOR 


Subgrade moisture behavior was definitely affected by highway design. The effects 
of improved shoulder width, drainage conditions, highway profile, traffic, and construc- 
tion procedure are given in this section. Instances of behavior and detailed case his- 
tories substantiating these findings are given elsewhere (6, 7, 12, 13). 


Improved Shoulder Width 


At all research sites observed, increasing shoulder width reduced the effects of 
runoff infiltration. Furthermore, in impervious sections with wide improved shoulders, 
the wet-dry interface or distinct transition point from fluctuation to more nearly stable 
subgrade moisture content occurred underneath the shoulders, whereas for pavements 
with open shoulders or narrow improved shoulders it occurred under the pavement, 
Less severe differential movements from more uniform moisture content were thus 
expected and observed for impervious pavement with wide shoulders. The wet-dry 
interface was found to occur approximately 5 to 7 ft from the outside edge of a covered 
area. Therefore, to keep the wet-dry interface from underneath the pavement, a 
shoulder width of at least 8 ft is required. 


Highway Drainage Conditions 


Drainage conditions or rainfall infiltration tendencies or both were found to be 
closely related to both shoulder slope and ditch design. Steeper shoulder slopes re- 
duced infiltration into the subgrade, and the elimination of ponding by the quick removal 
of surface runoff in ditches further reduced infiltration at most of the research sites 
where these drainage conditions existed. Gently sloping shoulder slopes were found to 
cause greater infiltration of runoff. However, ponding and improper ditch drainage, 
e.g., where ditches remained full for several days after a rain, caused measurable in- 
filtration despite good shoulder slope conditions. 

Nonuniform infiltration, producing differential moisture contents across the pave- 
ment section, was found to occur at many research sites. For undivided highways, this 
behavior usually occurred when the highways were constructed normal to an existing 
slope, with good drainage on 1 side and a tendency for ponding on the other. However, 
differential subgrade moisture contents from different infiltration rates were found 
under almost all 4-lane divided highways. It appears much more difficult to obtain 
good median drainage, and the median side of the divided highway usually experienced 
both higher maximum and average moisture contents and greater fluctuations at upper 
subgrade levels. 

Good drainage produced lower average subgrade moisture contents at the cost of 
higher variation and distinct wet-dry interface behavior. Poor drainage produced 
higher average moisture contents but less variation. 


Highway Profile 


Type of highway profile also affected subgrade moisture conditions; pavements on 
grade or in slight cuts usually exhibit both higher pavement ratings and more nearly 
uniform subgrade moisture conditions, particularly in the upper subgrade levels. Fill 
and transition sections usually had lower pavement and shoulder ratings than other sec- 
tions and also exhibited more moisture content variation, particularly at the pavement 
edges. 


Traffic Volume 


Despite several attempts during the study, significant correlations between traffic 
volume and subgrade moisture-related highway performance could not be established. 
Even when traffic volume was reduced to the general classifications of light, medium, 
and heavy, the only conclusion found was the obvious one that, after initial cracking or 
rutting or both of these rendered the pavement system pervious, higher traffic volumes 
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produced more rapid pavement deterioration. These data, coupled with observations 
that pervious bases were not usually found with degrees of saturation high enough to 
negate conventional wheel load distribution assumptions, form the basis for the con- 
clusion that subgrade moisture conditions are more important than traffic volumes in 
determining initial (and thus, to some extent, final) pavement performance. This con- 
clusion is also based on the assumption that at least a semirational method of con- 
sidering traffic loadings and pavement stresses was used in pavement system design. 


Current Oklahoma Highway Construction Procedure 


Most new construction in Oklahoma is done by stage methods under a long-term 
planning program. Contracts for grading and drainage structures are usually let and 
completed, and contracts for base and surfacing are let at some later time. This con- 
struction practice results in bringing the subgrade to its approximate final level and 
then in leaving it for an extended period before a covering or wearing surface is ap- 
plied. During the interim, several things occur. First, moisture contents in the pre- 
pared subgrade reach some equilibrium condition compatible with moisture contents in 
the natural subgrade below and with existing climatological conditions. This equilib- 
rium is likely to undergo seasonal changes, becoming wetter in general during winter 
months and drier during summer months, Also, precipitation in Oklahoma is usually 
infrequent but heavy when it does occur. As a result, considerable runoff occurs across 
the prepared subgrade and produces erosion damage. 

Another fairly common practice in Oklahoma is the use of prepared subgrades by 
farmers for movement of agricultural equipment either prohibited from or hazardous 
to on-pavement travel. If done for an extended period of time, this practice results in 
raveling of the subgrade during drier periods and rutting of the subgrade during wetter 
periods. 

As a result of these factors, the original subgrade level is changed and no longer 
remains suitable for use as a working surface. The base-surfacing contractor is, 
therefore, usually required to scarify and recompact at least the upper portions of sub- 
grade along the highway profile. Also, highway construction in Oklahoma is usually 
performed in warm spring and hot summer months. At this time the moisture in the 
uncovered subgrade is likely to be less to a considerable depth than the original com- 
paction moisture. Compaction of cohesive subgrades in Oklahoma is usually done dry 
of optimum because optimum compaction moisture is either slightly below or at the 
plastic limit of the material. The end result is to produce a subgrade with an average 
moisture content in the upper portions of at least 5 percent moisture content below the 
subgrade plastic limit. Base and impervious (at least initially) surfacing are immedi- 
ately applied over this material, and then subgrade moisture accumulation begins. 

The author is not criticizing the idea of staged construction because many advan- 
tages exist for its use, particularly in Oklahoma. However, other methods of staged 
construction should be considered that would provide better initial subgrade moisture 
conditions. Current construction practice appears to be aggravating the subgrade 
moisture problem in Oklahoma instead of minimizing it. 


EFFECTIVENESS OF HIGHWAY COMPONENTS IN PREVENTING AND 
RESISTING SUBGRADE MOISTURE CHANGES 


This section discusses various pavement system components currently being used 
in Oklahoma highway construction and the resistance of components to subgrade mois- 
ture variations. Recommendations for use of certain pavement components in particu- 
lar situations are also given. Additional descriptions of behavior, collected data, and 
case histories are given elsewhere (8, 9, 10, 11). 


General Philosophy for Highway Design on Expansive Oklahoma Subgrades 


Most published criteria for pavement design to resist effects of subgrade moisture 
are concerned with keeping moisture out of the subgrade. However, particular environ- 
mental conditions existing in the more populous areas of Oklahoma, coupled with 
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methods currently favored for all new construction, make it unlikely that construction 
moisture contents can be maintained for even short periods. 

Therefore, it appears logical to develop highway designs that allow subgrade mois- 
ture contents to increase to their equilibrium condition as quickly as possible and to 
stay there. The system should be designed to remain impervious after differential 
vertical and lateral expansion associated with Oklahoma subgrade moisture accumula- 
tion has occurred. Deterioration of Oklahoma highways must be prevented by not al- 
lowing the infiltration and evaporation cycle through pervious pavement systems to 
begin. 


Performance of Surfacing 


Current Oklahoma highway design procedures involve use of both flexible and rigid 
pavements. Type of surfacing or surface course used was found to have little, if any, 
initial effect on observed subgrade moisture behavior. Of more importance was 
whether the surface was pervious or impervious because the wearing surface itself 
simply serves to keep moisture out of the subgrade or else let it infiltrate through 
cracks, 

Rigid pavements were found to be more sensitive to longitudinal cracking than their 
flexible counterparts and were extremely sensitive to moisture infiltration through 
joints opened from thermal contraction and lateral subgrade expansion. However, rigid 
pavements were found to perform adequately for extended periods, after initial cracking 
and resulting infiltration had begun, and to maintain relatively good riding character - 
istics at subgrade moisture contents higher than those that were observed for flexible 
systems. The phenomena of pumping, widespread and severe cracking, subgrade shift- 
ing, and rapid deterioration did not usually occur until subgrade layers immediately 
under the pavement reached moisture contents considerably above the plastic limit and 
sometimes near the liquid limit. 

Flexible pavements or pavement systems consisted either of asphaltic concrete sur- 
facing over some other type of base or subbase (or both) or else asphaltic concrete sur- 
facing over some other type of asphaltic base or subbase. In general, better perfor- 
mance was obtained from multilayer asphaltic systems, as opposed to asphalt over 
nonasphalt systems. Flexible pavements were found to be highly resistant to cracking 
caused by small differential vertical movements and lateral subgrade expansion, es- 
pecially when more than 1 component of the system was composed of asphaltic ma- 
terials. However, flexible pavements were extremely sensitive to failure by loss of 
subgrade support when the moisture content of the upper subgrade material approached 
and exceeded the plastic limit. Failure was characterized by rutting, which was fol- 
lowed by pavement cracking. Cracks usually extended down to the subgrade material, 
and infiltration and evaporation of rainfall produced further deterioration and general 
cracking. Better observed performance of all-asphaltic systems probably results more 
from a thicker layer of asphaltic material that must be cracked to allow infiltration and 
evaporation than from any outstanding load-bearing or distributing characteristics of 
the material. 


Recommendations Concerning Type of Surfacing 


Portland cement concrete and asphaltic concrete surfacing are the only 2 choices 
available to the highway engineer. When pavements are constructed on expansive sub- 
grades, it is definitely recommended that the entire pavement system design be based 
on type of subgrade and other existing conditions rather than on use of a standard 
section. 

Initial resistance to cracking from vertical and lateral subgrade expansion was 
lower in rigid pavements than in flexible pavements. However, rigid pavements were 
found to perform adequately (even though cracked) at subgrade moisture contents well 
above the plastic limit because of their ability to transmit wheel loadings over a wide 
area. The tendency of rigid pavements to crack with subgrade volume change is di- 
rectly responsible for increasing moisture contents; and, once surface cracking has 
occurred, the deterioration of rigid pavement is only a matter of time. Therefore, 
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rigid pavements should not be used unless underlaid by a nonexpansive, flexible, im- 
pervious material. Load-carrying capacities of rigid pavement should be evaluated by 
strength tests, assuming that there is a subgrade moisture content of 1.1 to 1.3 times 
its plastic limit, The ability of rigid pavement to effectively carry traffic at subgrade 
moisture contents above the plastic limit should indicate its use in areas where poor 
drainage is encountered because rigid pavement normally produces higher subgrade 
moisture conditions. A correctly designed rigid pavement with improved shoulders 
over a sand-asphalt layer is probably the best design for general use on expansive 
subgrades, 

Initial cracking from vertical and lateral subgrade expansion was resisted better by 
flexible pavement than by rigid pavement. However, when moisture contents approached 
the equilibrium value of 1.1 to 1.3 times the subgrade plastic limit, heavy traffic pro- 
duced rutting and initial pavement cracking. The inability of observed flexible pave- 
ment systems to carry heavy traffic loadings at subgrade moisture contents near the 
plastic limit is thought not to be an indictment of this design but a reflection of the in- 
ability of current Oklahoma Department of Highways desien techniques to predict sub- 
grade strength at these moisture contents. Incorporation of revised soil-testing pro- 
cedures to adequately determine subgrade strengths at moisture contents 1.1 to 1.3 
times the subgrade plastic limit and to consider volume change behavior should defi- 
nitely improve the chances for satisfactory long-term flexible pavement performance. 
Also, flexible pavement construction on expansive subgrades should definitely include 
the use of a nonexpansive, flexible, and impervious base, subbase, or membrane 
component, 

Subgrade strength appears to be more important in flexible pavement design, and, 
therefore, these systems should give better performance where good drainage condi- 
tions are provided because they are also less susceptible to cracking from volume 
change produced by good drainage and would benefit materially from lower average 
moisture contents and resulting higher subgrade strengths. 


Performance of Base Courses 


Six types of base material were encountered at field research sites that are com- 
monly used in Oklahoma base construction: sand cushion, hot sand asphalt, asphaltic 
black base, soil cement, stabilized aggregate, and select material. Exact definitions 
of these common highway materials are available elsewhere (8, 9). Sand cushion and 
select material were encountered as base courses for rigid pavement. 

Because of observed behavior, the use of sand cushions has been deleted from 
Oklahoma highway design criteria. Sand cushions were found to act as water reser- 
voirs and distribution systems, catching water that infiltrated from the shoulders and 
through rigid pavement cracks and joints and feeding it uniformly over the subgrade. 
As a result, rigid pavements on sand cushions usually experienced only small differ- 
ential vertical movements. However, continued feeding of water to the subgrade re- 
sulted in lateral subgrade expansion and longitudinal pavement cracking plus upper sub- 
grade moisture contents being at or near the subgrade liquid limit. 

Hot sand asphalt was used as a base course at many highway research sites and was 
found to form an excellent impervious layer. The plastic properties of asphalt that 
were coupled with the fineness of the mineral particles allowed this material to resist 
the effects of lateral subgrade expansion and vertical differential movement without 
cracking and becoming pervious, at least better than any other type of material previ- 
ously encountered. 

Several research sites contained asphaltic black base underneath asphaltic concrete. 
Observed performance indicates that this material is probably the second-best base 
type tested because it is impervious and possesses good flexibility. However, its per- 
formance is not thought to be as good as that of sand asphalt because the larger size of 
aggregate particles used probably gives this material less flexibility and also presents 
the possibility of larger interconnected voids. Nevertheless, asphaltic black base was 
found to perform well at the majority of sites where it was encountered. 
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The soil-cement base courses that were encountered were used primarily in the 
construction of improved shoulders. Soil cement is relatively rigid compared to more 
flexible asphaltic materials and may suffer some initial shrinkage cracking. These 
cracks expand when moisture accumulation and resulting vertical and lateral subgrade 
expansion occur. Cracking is almost immediately reflected through the thinly surfaced 
shoulder, and this random cracking may be observed on almost all Oklahoma highway 
shoulders applied over soil-cement bases. Once the surfacing has cracked, water will 
enter and infiltrate through the soil cement into the subgrade, and additional volume 
change will occur and consequently open the joint between pavement and shoulder. The 
use of soil cement directly on any subgrade of even suspected expansiveness does not 
appear to be advantageous. 

Stabilized aggregate is a mechanically stabilized material consisting of blended 
course aggregate, sand, mineral filler, and soil binder. The resulting product, although 
highly variable and depending on locally available materials, is nevertheless intended 
to provide a well-graded and densely compacted layer with reasonable strength prop- 
erties. Stabilized aggregate base courses were found to perform satisfactorily at sites 
where little lateral subgrade expansion was noted to occur. At these sites, they pro- 
vided (by virtue of their density and fine content) a relatively impervious barrier to in- 
filtration, both from the shoulders and through the pavement section, However, sta- 
bilized aggregate possesses little tensile strength; and, at sites where appreciable 
lateral expansion was noted, this base did not stop infiltration and evaporation of water 
through pervious pavement surfaces. 

Select material is similar to stabilized aggregate but is usually considered to be 
naturally occurring and is governed by somewhat different specifications. On the 
whole, select material was observed to behave not quite so well as stabilized aggregate 
but better than sand cushions. 

Most of the correlations concerning base material performance are available else- 
where (8), but it should be noted that, comparing subgrade moisture contents under 
pavement centerlines for research site locations on more or less uniform A-6 or A-7 
subgrade and on pavement rated as being excellent or good, 18 of 23 sites with pervious 
bases had higher moisture contents in upper subgrade levels than in lower levels. On 
the other hand, 7 of 11 similar pavements with impervious bases had lower subgrade 
moisture contents in the upper subgrade under their centerlines than at lower subgrade 
levels. Evaluation of amount of subgrade moisture variation under pervious and im- 
pervious bases, although subjective to some degree, nevertheless indicated that less 
variation with time existed in upper subgrade levels under impervious bases and that 
relative magnitude of variations was also lower. 

The subgrade moisture contents for almost all of these sections were in the range 
of 1.1 to 1.3 times the plastic limit, with the upper subgrade moisture levels under the 
7 impervious base sites being closer to 1.1 than to 1.3 times the plastic limit. On the 
other hand, average moisture content-plastic limit ratios for sites on pervious bases 
approached (and sometimes exceeded) the 1.3 value. These data, plus other informa- 
tion (8), indicate that a reasonable amount of moisture infiltrates into the subgrade 


from pervious base courses, even when the pavement is still in acceptable condition. 


Recommendations Concerning Highway Base Courses 


Base courses used in Oklahoma highway construction on expansive subgrades 
should, if possible, be nonexpansive, flexible, and impervious. Asphaltic types of base 
materials appear to provide these qualities better than any other materials encountered. 
If use of other base material is contemplated, the subgrade should be protected by use 
of a flexible, impervious subbase or else a flexible, impervious membrane located 
somewhere in the pavement system. If an impervious layer is provided in some other 
portion of the pavement system, mechanically stabilized aggregate or even select ma- 
terial may be used, as long as it does not become highly saturated. Soil cement should 
definitely not be used as a base material unless it is protected from lateral subgrade 
expansion by an intervening nonexpansive, flexible, and impervious layer. 
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Recommendations Concerning Highway Subbase Courses 


Subbase components are used primarily as a means of distributing loads to the sub- 
grade, but they can also act as barriers against water infiltration where pervious base 
courses are used. Oklahoma subbases are usually constructed of select material, 
chemically treated layers of the subgrade, or sand asphalt. 

Select material subbases, when used in flexible pavement construction, obviously 
reduce the total required thickness of asphaltic layers. Even though economics dictates 
the use of locally available materials whenever possible, especially if they do not re- 
quire treatment, it should be remembered that select material has minimal tensile 
strength and is not completely impervious. Thus, select material subbases should not 
be used unless some other portion of the pavement system (besides the wearing sur- 
face) contains desired nonexpansive, impervious, and flexible properties. However, the 
dense and relatively impervious qualities of select subbase are helpful in reducing ef- 
fects of infiltration from outside the pavement system and also in protecting the base 
and wearing surface components from subgrade moisture and volume changes by pro- 
viding a cushioning effect, 

Lime-treated layers of cohesive subgrade have also been used as subbases, normally 
with lime contents corresponding to modification optimum for the material. The lime 
contents used for treatment normally result in reduction of plasticity and give an in- 
crease in soil workability, The result is to produce a treated material with properties 
similar to that of the select material described previously. 

Sand asphalt possesses necessary imperviousness and flexibility if these qualities 
are desired in a subbase, and its use should allow inclusion of pervious base materials 
in highway design. Soil cement should not be used as a subbase if it is expected to re- 
main reasonably impervious, and, if used, severe cracking should be expected. 


Recommendations Concerning Improved Shoulders 


Improved shoulders were found to provide a method for reducing infiltration of sur- 
face runoff into the base course, subbase, or subgrade or all three, particularly at 
pavement edges, and thus for producing more nearly uniform moisture variations over 
the width of the pavement. However, definite improvement in subgrade moisture condi- 
tions underneath pavement could be obtained by several modifications in shoulder de- 
sign and construction technique. Because of lower traffic loadings carried by improved 
shoulders and also the staged construction common in Oklahoma, shoulders are often 
built after the pavement section has been constructed and have thinner base layers that 
are often of a different material. The net effect of this construction technique is to 
place a vertical plane of weakness between pavement and shoulder. Lateral and differ- 
ential vertical expansion of cohesive subgrade tends to cause separation along this plane 
of weakness, provides a channel for entrance of surface runoff, and produces larger 
moisture variations under the pavement edge. Also, the wet-dry interface that forms 
under the shoulder causes shrinking and swelling of subgrade material, which develops 
flexural stresses along this plane of weakness. 

Improved shoulders at least 8 ft in width should be constructed for all pavements on 
expansive subgrades. Base or subbase under the shoulders, or both of these, should be 
continuous with that of the pavement section and should be applied in 1 continuous 
shoulder-to-shoulder lift. Recommendations concerning type of base and subbase ma- 
terials under shoulders are the same as for wearing surfaces. 

Whatever the desired design shoulder width, shoulders on each side of the pavement 
section should be of equal width. If both shoulders are of equal width and the width is 
at least equal to 8 ft, a more nearly symmetrical subgrade moisture profile will result 
underneath the pavement and more nearly uniform moisture increases and decreases 
will occur. 


Recommendations Concerning Highway Drainage Conditions 


The varied topography encountered along a highway profile makes specific drainage 
recommendations impracticable. However, it is extremely important that existing 
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drainage conditions be the same on each side of a particular highway section to produce 
more nearly symmetrical moisture profiles, although conditions themselves may change 
along the profile. Special care should be given in transition sections and sections cut 
through slopes to provide equal drainage conditions on either side of the pavement. 
Also, despite numerous problems likely to be encountered, drainage conditions for 
4-lane divided highways should be designed so that both the median and the outside 
shoulders drain, As a general rule, flexible pavements appear to be more suitable for 
use under good drainage conditions, whereas rigid pavements appear to be better 

suited for poor drainage conditions. For intermediate conditions, correctly designed 
pavements of either type should work satisfactorily. 


Recommendations Concerning Highway Profile 


The diverse topography encountered along highway profiles makes specific recom- 
mendations impractical. However, research findings definitely indicate that better 
highway performance is obtained in sections on grade and in slight cuts. Where pos- 
sible, the highway profile on expansive subgrades should be held to these sections. Any 
fill or transition sections or both should be partially constructed of select material or 
lime-modified layers of the subgrade, such that a 5- to 7-ft depth below the bottom of 
the pavement system is essentially nonexpansive. An alternate procedure might be to 
compact the fill as closely as possible to natural moisture contents of existing soils 
and cover the entire fill with an asphaltic membrane. Whenever possible, high fills 
should be avoided, and, in any case, despite the previously mentioned treatments, 
shoulder settlements should be expected. 


Recommendations Concerning Current Staged-Construction Procedures 


If current staged-construction procedures (as described previously) are to be used, 
several slight modifications will improve moisture conditions in the subgrade. One 
method might be to compact the subgrade to a level above design final grade. This 
additional material will reduce subgrade drying during hotter seasons and give natural 
moisture contents more nearly approximating compaction specifications. In addition, 
this added layer will help to compensate for thickness lost through erosion and settle- 
ment, The prepared subgrade should be adequately barricaded so that agricultural 
traffic and other traffic may be kept off during the interim between subgrade finishing 
and base-surfacing application, The base-surfacing contractor could then cut the sub- 
grade to final grade and immediately apply base and surfacing. If this procedure is not 
feasible, consideration should be given to provide extra rolling time (and thus money) 
to compact upper subgrade layers wet of optimum. However, required densities may 
be difficult to obtain as compaction under these conditions will normally occur at or 
above the plastic limit. Nevertheless, higher initial moisture contents and a clay 
particle orientation less conducive to volume change would result and may be worth 
the additional effort. 

Another alternative procedure would be to cover the prepared subgrade with an 
asphaltic membrane and allow it to reach equilibrium moisture conditions over a a- 
year period or longer between subgrade completion and base-surfacing application. 
Still another alternative procedure would be the use of the deep-plow lime treatment of 
upper subgrade layers to produce a less expansive buffer layer, which would maintain 
more nearly constant moisture conditions in lower portions of the subgrade. This pro- 
cedure has been attempted experimentally by the Research and Development Division, 
Oklahoma Department of Highways (15). 


Recommendations Concerning Revised Staged-Construction Procedures 


Some thought should also be given to a more radical revision of the staged- 
construction process to ensure better long-term subgrade moisture conditions. The 
philosophy and methods by which current Oklahoma highways are constructed dictate 
that the entire base and surfacing courses be applied at one time. However, the ulti- 
mate traffic loading for the pavement system is not usually encountered until late in its 
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design life. Initially, traffic is likely to be substantially below expected maximum 
values, and the section is perhaps overdesigned. The possibility of traffic considered 
in design being greater than that during initial years after construction and establish- 
ment of subgrade moisture equilibrium suggests that alternative methods of staged 
construction be used. 

One method would be to divide construction into 3 stages. Initially, the subgrade 
would be prepared in a manner recommended previously. Then, after the normal time 
lag, base material would be applied, but only part of the total surfacing thickness would 
be applied. Surfacing applied at this time should be adequate to carry current traffic 
but should be less than the ultimate surfacing thickness recommended for the section. 
Within a reasonable length of time after initial surfacing is applied and the section is 
opened to traffic, moisture equilibrium should occur underneath the pavement system, 
and its subgrade will be just beginning to start the infiltration and evaporation cycle 
through the now pervious pavement structure. At this time, usually somewhere between 
2 and 5 years after initial construction, the final surfacing course should be applied to 
bring the pavement section up to design thickness and effectively seal the surface from 
infiltration and evaporation. Moisture equilibrium should have been achieved by this 
time, and, as the accumulation phase is completed, sections should be produced where 
adverse moisture variations would be minimized. If subgrade soil conditions at mois- 
ture equilibrium values have been correctly anticipated, excellent highway performance 
should be achieved, and future maintenance costs should be markedly reduced. During 
the design life of the highway, at least 1 less overlay should be needed, and this in itself 
would reduce total highway cost. 

An alternate procedure would be to use current staged-construction methods but 
then to apply surfacing in 2 stages. This procedure should cause more initial moisture 
accumulation and relative movement and might require the final surfacing course ear- 
lier than the procedure mentioned previously. However, it should also increase sub- 
grade moisture resistance of the final section. 

Still a third alternate might be to extend the initial staged-construction phase to in- 
clude base and temporary surfacing, with intermediate or final surfacing courses or 
both courses to be applied at later dates. This technique might be applicable when 
turnkey or non-staged-construction projects are anticipated because of considerations 
dictating rapid availability of the highway to traffic. The sections could be completed 
quickly, but provision would be available to counteract adverse subgrade moisture con- 
ditions, which would of necessity be built into a turnkey project. 


SUMMARY 


This paper has described subgrade moisture conditions in Oklahoma as determined 
by a 6-year research study. Effectiveness of various highway components in preventing 
and resisting adverse subgrade moisture conditions was also discussed. Recommenda- 
tions concerning the need for more comprehensive subgrade soil-testing procedures, 
for use of certain pavement components in particular situations, and for revision of cur- 
rent construction practice to enhance subgrade moisture resistance of pavement sys- 
tems have also been described. It is suggested that, whenever possible, recommenda- 


tions made here be applied in routine design and construction of highways on expansive 
Oklahoma soils. 
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PREFABRICATED SUBSURFACE DRAINS 


Kent A. Healy and Richard P. Long, Department of Civil Engineering, 
University of Connecticut 


The design of a prefabricated subsurface drain system based on well screen 
criteria is presented. This drain is fabricated by using synthetic cloth and 
a channelized core that ensures proper filtration and permeability. Labo- 
ratory and full-scale field test results are given and showthe prefabricated 
subsurface drains to be an effective substitute for conventional mineral 
aggregate subsurface drains. Mineral aggregate filters are reviewed. The 
difficulties inherent in their design and construction and the advantages of 
prefabricated subsurface drains are discussed. 


eSUBSURFACE DRAINS have been used for many years to remove excess water from 
the ground to improve crop growth, strengthen pavement foundations, stabilize slopes, 
and reduce water pressure against retaining walls. The water is often carried away 
from a site by a perforated pipe, and a filter material must be provided between the 
soil and the pipe that will retain the soil while allowing the free flow of water into the 
pipe. The conventional subsurface drain uses a filter of mineral aggregate that has 
been graded to match the soil (1). This type of filter performs satisfactorily only if 
designed and constructed carefully. 

A subsurface drain system that uses synthetic materials and that fulfills the re- 
quirements of filtration and water flow has been developed. Because this system can 
be prefabricated in lengths that are easily handled and installed in the field, many of 
the construction problems of mineral aggregate drains are eliminated. 


DESCRIPTION OF PREFABRICATED SUBSURFACE DRAINS 


The prefabricated drain is shown in Figure 1 and consists of a slotted pipe, a chan- 
nelized vertical core inserted into the pipe slot, and a fine-mesh filter cloth enclosing 
the pipe and the core. The cloth retains the soil and keeps the core channels open, 
The groundwater drains through the cloth, down the channels, into the pipe, and away 
from the site, 

The selection of the filter cloth is based on criteria established in the 1930's for 
well screens (2). A properly designed well screen retains only the coarsest particles 
of a soil; these particles will retain the finer particles. 

The openings in the screen should be just small enough to retain the coarsest 20 
percent of the particles (D,,) in a well-graded soil and the coarsest 40 percent (D,,) in 
a uniform soil (3). The percentage of open area of the screen should be approximately 
the same as the porosity of the soil to prevent restriction of water flow. Well screens 
meeting these criteria have remained effective for many years in soils ranging from 
fine sand to gravel. 

The filter cloth used on the prefabricated drains functions in the same way as a well 
screen. An advantage of using cloth as a filter is its high permeability even with small 
openings. Cloth with mesh openings between 0.075 and 0.150 mm and a minimum 15 
percent open area will retain most soils that can be effectively drained by gravity and 
will not restrict the flow of water from the soil. 


Sponsored by Committee on Subsurface Drainage and presented at the 50th Annual Meeting. 
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Materials Used 
The materials used in the prefabricated surface drains are given as follows: 


1, Filter cloth—Two types of cloth have been tested and used as filters—a nylon 
chiffon and a polyester butterfly (Fig. 2). Both materials have good tensile strength 
and resist decomposition in the ground. The chiffon has a mesh opening of 0.15 mm 
and a 45 percent open area, The butterfly has a mesh opening of 0.075 mm and a 15 
percent open area. Long-term model tests have shown that both cloths will filter soils 
ranging from a fine silty sand to a glacial till with no sign of soil erosion or clogging. 
The butterfly cloth with 15 percent open area has adequate permeability for soils con- 
taining silt. The chiffon having a greater percent of open area is preferable for drain- 
ing more permeable soil. Both cloths are more permeable than a 2-in. layer of coarse 
gravel. The gradation of particle sizes and permeability of soils successfully filtered 
by nylon chiffon cloth are shown in Figure 3. These soils, ranging from a medium sand 
to a glacial till with a high percentage of silt, were tested in the laboratory under a 
hydraulic gradient much higher than that occurring in the field. 

2. Core—The core must support the cloth and provide channels large enough to 
carry the water into the pipe as quickly as it flows from the soil. The following mate- 
rials have been used as cores in field-test installations and are shown in Figure 4: 
Type 1, expanded aluminum sheet, purchased from U.S. Steel, having the commercial 
name Armorweave; and Type 2, vinyl tube fencing, purchased from Sears, Roebuck and 
Co, The core materials were tested for crushing strength and water-carrying ability. 
The expanded aluminum showed little deformation under a pressure of 200 lb/sq in., 
and the vinyl tube fencing had a crushing strength of 160 Ib/sq in. Intrusion of the cloth 
into the channels decreases the area through which the water can flow and is a possible 
source of malfunction. Tests with the vinyl 
tube fencing showed that, under an earth 
pressure of 4.0 kip/sq ft, the cloth intruded 
less than 0.025 in. and that a 1-ft wide 
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Figure 1. Basic design of prefabricated subsurface Figure 2. Weave characteristics and porosity of cloth 
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Figure 3. Soils successfully filtered by nylon chiffon. 


section of the drain, when vertical, could carry 2.5 gal/min. This is sufficient to drain 
a soil with a permeability of 5 x 10-?ft/min. Tests on the expanded aluminum core 
showed similar flow characteristics under stress. 

3. Pipe—A hard plastic pipe with a 4-in, diameter and a 1/8-in. wall thickness has 


Type l-Expanded Aluminum 


BINDING 


Figure 4. Core configuration. 


been used for the field sections. The 
slot was cut on a table saw. This 
pipe is available in 10-ft lengths from 
building supply companies and has a 
crushing strength of 0.9 kip/ft. 


Installation Methods 


Prefabricated drains are easily 
installed on slopes or in trenches, as 
shown in Figure 5. Drain sections 
10 ft long weigh less than 20 lb and 
can be placed and connected with slip 
couplings from the ground surface, 
Trenches need be only wide enough to 
receive the drain, thereby making 
attractive the use of trenching ma- 
chines in suitable soil. 

The core is flexible enough to be 
pressed tightly against 1 face of the 
exposed natural soil, as shown in 
Figure 5, In this manner the more 
permeable strata in the native soil can 
be drained quickly. The type of back- 
fill depends on the drainage desired 
at the site. If most of the water 
enters the drain from 1 side, as ona 
slope, the native soil may be used as 
backfill against the opposite face. 
When the drain receives water from 
both sides, sand may be used as back- 
fill, This ensures that the natural 
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drainage channels in the soil are not blocked 
and makes the compaction easier in atrench 
section. The short flow path through the 
sand to the core of the drain will not im- 
pede the free movement of water. 

The prefabrication of the drain sections 
ensures proper operation even when in- 
stalled by people who are unfamiliar with 
filter principles. The drain sections can 
be fabricated in any height and length to 


NATIVE SOIL 


CORE 


suit the installation. SAND 
BACKFILL 
CHUTED 
FIELD TESTS INTO 
, TRENCH 
Sites 


Prefabricated drains were installed in 
2 small wet areas in 1968 and performed 
well. In the spring of 1969, 600 ft of drain 
were installed to stabilize a cut slope that AY. 4" PIPE 
was sloughing because of excess water. 
This installation is described in detail in wee 
this paper. 

During the late summer and early fall 
of 1970, installations were made at 3 other 
sites. These 3 installations are mentioned 
especially to illustrate other applications for the prefabricated underdrains. Field 
measurements on them are continuing. One installation surrounds a septic system 
leach field to control the groundwater in the vicinity. Another installation was placed 
to intercept the groundwater flowing into a lot on which a home is to be built. The third 
installation was placed beside a road to control frost heave by lowering the local 
groundwater table. 


Figure 5. Prefabricated drain installed in trench. 


Description of Unstable Slope 


The slope, a plan of which is shown in Figure 6, is the northwest side of a drumlin, 
located on the University of Connecticut campus and formed when the hill was cut back 
from a natural slope of 1 on 3.3 (17 deg) to a 1 on 2 slope to allow the placement of a 
sanitary line and sidewalk. The slope started sloughing after the first heavy rain, and 
it was a continual maintenance problem to keep the walk clear of mud in the spring and 
ice in the winter. The soil in its natural state is a dense, well-graded glacial till with 
particles varying from cobbles to clay size. Disturbed samples of the soil have a per- 
meability of approximately 1 x 10-° ft/min, as measured by a falling head permeameter. 
Slow direct shear tests showed an effective stress friction angle of 41 deg. The natural 
undisturbed soil is slightly cemented and contains numerous small channels parallel to 
the surface that seep water below the water table in an open cut. 


Installation 


In July 1969, 2 lines of prefabricated drains with type 1 core and butterfly cloth 
were installed along the slope, as shown in Figure 6, It had originally been planned to 
install the upper drain line by cutting a berm with a bulldozer. However, at the time of 
installation the slope was too wet, and a berm had to be dug by backhoe from below. 
Figure 7 shows a typical cross section and the method used to place the drains. A 
trench was dug for the lower drain and was backfilled with sand for ease of compaction. 
In the very wet areas, the upper drain sections were partially backfilled with the back- 
hoe immediately after placement to prevent local sliding. Figure 8 shows the upper 
drain partially backfilled. Final backfilling and grading were completed by a bulldozer, 
as shown in Figure 9. The lower trench installation is shown in Figure 10. 
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SIDEWALK 


Figure 6. Plan of field installation. 


In the fall of 1969, heavy rains caused surface erosion, and in 2 areas, marked A and 
B in Figure 6, the natural drainage channels were such that water was exiting under the 
upper drain, causing sloughing below. An additional 20-ft length of drain was installed 
by hand at the south end of the area B halfway up the slope, and this portion was 


Upper Drain Installation Sequence 
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|. BERM CUT WITH 2. DRAIN IN PLACE 


3. BACKFILLING WITH 
BACKHOE NATIVE SOIL 


Cross Section of Slope 


+ 
UPPER DRAIN 


SIDE WALK 
4! 
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Figure 7. Cross section of field installation. 
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stabilized. In the spring of 1970, 90 ft of 
drain with type 2 core and chiffon filter 
cloth was installed with a backhoe and 
bulldozer in the remaining portion of 
area B, and 20 ft of drain was installed in 
area A. Observation pipes were installed 
at points marked 1, 2, and 3 in Figure 6, 
and the water flow in the drain pipes can 
be measured at these points by using 
calibrated probes. 


Evaluation of Field Installation 


Some piezometers were installed, but 
they did not reflect apparent water condi- Figure 8. Upper drain partially backfilled. 
tions accurately. This may be due to the 
nonhomogeneous permeability character- 
istics of the soil, The evaluation of the field installation has therefore, been based on 
the overall stability of the slope and measurements of water flow out of the drains. 

The slope, which was unstable over essentially its whole length, has been stabilized. 
Surface seeping has been almost eliminated, and the slope surface drys up within a few 
days after a heavy rain. 

Water flow from the drains has been monitored continually since installation. The 
250-ft upper drain removes water from the soil at a rate of 2 gal/min during wet periods. 
Calculations, assuming a hydraulic gradient in the soil of 0.2, show the field perme- 
ability to be 1.6 x 10-3 ft/min. This increase of 1,000 times over the permeability from 
lab tests on disturbed samples may be due to the natural channels occurring in the de- 
posit. There has been no indication of fines being removed from the soil by the drain, 
indicating that the cloth is filtering properly. 


Costs 


The upper drain was installed where it would have been nearly impossible to install 
a conventional trench drain, so that cost comparisons are difficult to make, but the 
prefabricated drains appear to be competitive with conventional drain systems. The 
material cost of the prefabricated drain using vinyl tube core and chiffon is less than 
$1.35/ft for a 4-ft deep section. The main drainage core can be provided in almost any 
height and length at a cost of approximately $0.25/sq ft. 


COMPARISON OF MINERAL AGGREGATE AND 
PREFABRICATED SUBSURFACE DRAINS 


Most subsurface drains are constructed from mineral aggregates, and a comparison 
with prefabricated drains is appropriate. 


Figure 9. Dozer backfilling upper drain. Figure 10. Installing lower drain in trench. 
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Filter Requirements 


Tests run by Bertram in 1940 (4) and at the U.S. Army Engineer Waterways Experi- 
ment Station (5) resulted in the development of the following criteria for mineral aggre- 
gate filters: 


1, Preventing continuous movement of soil particles requires that the effective 
pore size (assumed to be 1/5 D,, of the filter) be smaller than the coarsest 15 percent 
(D,,) of the soil being drained. This is normally expressed as 


D,,; (filter) 


D,.teuy * 5 


Deo (filter) 


D.y (soil) £ & 


2. Preventing restriction of water flow by the filter requires 


D,; (filter) 


D,; (soil) . ee 


Preventing movement of the filter particles into the pipe, if a perforated pipe is used to 
remove the water from the filter, requires 


Dg; (filter) 


Ppipe-opening size 


In the prefabricated subsurface drain, all filtration is accomplished by the cloth with 
openings constituting at least 15 percent of the area and having a size between 0.075 
and 0.150 mm. 


Design and Construction 


In many Situations the criteria for mineral aggregate filters must be applied with 
great care. The Vicksburg criteria implicitly assume well-graded soil (1). If the soil 
to be filtered and drained is gap-graded, the number of large particles may be insuf- 
ficient to prevent movement of the smaller particles, and the filter must be designed to 
retain a size smaller than the coarsest 15 percent of the soil particles. 

Filters placed against soil deposits, whose gradation varies from point to point, 
must be designed to hold the finest particles in place, and a graded mineral aggregate 
filter may be required to allow free drainage. 

Some of the important points (1) in constructing the mineral aggregate filter are as 
follows: 


1. Filter materials must be handled and placed with care to avoid segregation and 
contamination; 

2. The filter must be well compacted to reduce the possibility of dropping fines of 
the filter through void spaces; and 

3. A single improperly constructed portion of the filter can lead to failure of the 
drainage system. 


Construction control for the prefabricated subsurface drain is less demanding. The 
prefabrication ensures that the system can be easily and correctly installed by person- 
nel unfamiliar with filter criteria. 

CONCLUSIONS 


Laboratory and field tests have indicated the following: 
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1. A fine mesh cloth is suitable as an effective filter for a wide range of soil types; 

2. A thin channelized core allows free movement of water into the outlet pipe; 

3. Prefabricated subsurface drains are easily handled and installed in the field and 
allow placement where conventional drains would be difficult to construct; and 

4, Prefabricated subsurface drains are economically competitive with conventional 
mineral aggregate systems. 
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GRAVITY FLOW TO EXCAVATIONS AND DRAINAGE 
TRENCHES IN LAYERED AQUIFERS 


Francis G. McLean, Westenhoff and Novick, Inc., Chicago; and 
Raymond J. Krizek, The Technological Institute, Northwestern University 


The finite-element method is used to analyze the problem of steady-state 
gravity flow to typical excavations and drainage trenches in layered aquifers. 
Dimensionless flow quantities and information on the location of the phreatic 
surface are presented as functions of the relative material permeabilities 
and the geometric configuration of the soil profile and excavation size or 
drainage trench position. These results are then applied to an example 
problem to select a satisfactory configuration of subsurface drains for an 
actual depressed highway profile. 


eALTHOUGH the topic of plane flow through layered systems has been investigated ex- 
tensively by workers in the field of agricultural drainage (12), a large portion of the 
effort has been directed toward the solution of confined flow problems (9, 12, 13, 14). 
Relatively little attention has been given to gravity flow systems (such as flow to ex- 
cavations and drainage trenches), which are commonly encountered in civil engineering 
practice in general and in highway engineering in particular. Accordingly, the multiple 
aquifer systems shown in Figure 1 were studied to gain insight regarding the interac- 
tion effects of the various layers when subjected to conditions of steady-state gravity 
flow. The single excavation without drains (Fig. 1a) is typical of a general system with 
wide application on construction sites or in highway cuts, whereas the configuration 
with drains (Fig. 1b), although broad in use, is limited for illustrative purposes to a 
particular combination of material permeabilities that are representative of those 
found in a portion of a large highway project. The latter more specialized situation is 
used to test the concepts developed during the more general study of the single excava- 
tion. Because the combination of a permeable boundary at the trench side walls, the 
free surface, and the layered materials (with the resulting complex configuration of the 
free surface) makes these problems very difficult and generally tedious to solve by 
ordinary methods, the finite-element method (4, 6, 17, 20, 21, 22) was chosen for use 
in this investigation. 


EVALUATION OF FINITE-ELEMENT METHOD 


The advantages and disadvantages, as well as the accuracy, of the finite-element 
method has been well documented by workers in the area of structural analysis. Some 
typical studies include the consideration of various formulations (5, 8) and the descrip- 
tion of various structural systems by different types of discrete elements (15, 21). 
Much has been published about the methods for obtaining solutions to the associated 
system of simultaneous equations, the conditioning thereof, and accuracy and error 
analysis of the procedures and answers. Results obtained by the finite-element 
method, as applied to the 'quasi-harmonic" problem (which includes seepage), have 
been shown (19, 20, 21) to be comparable with those obtained by finite-difference 
methods and closed-form solutions. Also, a comparison has been made (6, 17) be- 
tween results obtained by a finite-element analysis and those obtained by Casagrande 


Sponsored by Committee on Subsurface Drainage. 
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Figure 1. Typical geometric configurations considered. 


(2) for the location of the free surface and the exit point for steady-state flow through 
a dam that rests on an impermeable base. The effect of mesh size on the solution for 
a geometrically simple, confined-flow problem has also been reported elsewhere (11). 

The accuracy of the finite-element program (essentially that described by Taylor 
and Brown, 17) used in this study was ascertained by making several checks of flow 
quantities, seepage pressures, and free surface locations for various axisymmetric 
and plane flow cases. In the first case, an electric analog model of a trench, which 
partially penetrates a homogeneous, isotropic aquifer, was used (18) to check the flow 
quantity and potential distribution obtained by the finite-element method. The results 
shown in Figure 2 indicate very good agreement. In another comparison, 3 cases of a 
well, which partially penetrates an unconfined, homogeneous, isotropic aquifer under- 
laid by an impermeable stratum, were analyzed by the finite-element method. Results 
in all cases differed by less than 10 percent from those obtained (1) by relaxation and 
the methods of Kozeny and Forchheimer. ~ 

The problem of gravity flow through a dam on an impermeable base has been 
studied by relaxation (16), by the hodograph method (13), by an iteration scheme (7), 
and by use of a flow model (3). Hence, this problem affords an excellent opportunity 
for comparison of solutions. 
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Figure 3a shows the finite-element mesh that was constructed to describe this 
problem, and Table 1 gives the geometric information and the resulting flow and pres- 
sure characteristics for the specific cases considered. The maximum number of itera- 
tions allowed for each case was 15, and, if completed, they required about 200 sec of 
central processor time on a CDC 6400 computer. The specified beta quantity is used 
in the program as an under-relaxation factor for free-surface adjustments between 
iterations, and the maximum pressure variation (associated with particular values of 
beta and tolerance) for the nodes along the free surface is given in Table 1. Also shown 
are the dimensionless flow parameters calculated from the finite-element analysis and 
the Dupuit assumption, which, according to Muskat (13), yields values almost identical 
to those obtained by the hodograph method. The first trial and final locations for the 
free surface in each case are shown in Figure 3b, and Figure 3c shows the results for 
a problem whose geometry was chosen to compare with that of the experimental model 
investigated by Chapman (3). The solution obtained (10) by the method of finite differ- 
ences is also shown. Because of the excessive computer time required, the finite- 
element solution was not allowed to attain the best possible free-surface location, nor 

was a better first-trial free surface 

tried; however, all points that lie 

between the solution curves were 

aainsass moving upward when the solution 

* 420 elements SHG Sent bi Or Tit _ was terminated. The single point 
SS ee lying above the Chapman curve near 

“#i" the tail water was observed to be 

E = essentially stable during the last 2 
iterations. 

For the cases considered here, 
om ++ = the convergence of the solution de- 
- termined by the finite-element 

= _ idealization and the convergence of 
the free-surface iterative procedure 
are of concern. The latter may be 
evaluated for each specific solution 
(a) Finite Element Mesh for Flow Through a Vertical- faced Dam by observing the residual pressures 
calculated for the free-surface 
nodes, whereas the former may be 
assessed by comparing results ob- 
tained from several different ideal- 
izations for any given problem, In 
+ Hy general, as the size of the elements 


r L *| 


= Te ries sunleee.” : | becomes small and the number of 
—— Final Position of Free Surface | elements becomes large, conver- 

Pp AES 7 gence is ensured for most condi- 
te L=20 | tions (21). For the several hundred 


computer solutions performed during 
the course of this and associated 
studies, satisfactory convergence 

Sco duteemsecetnanioaen was obtained through the exercise 


(b) First Trial and Final Positions of Free Surface for Various Head Differences 


---- Jaeger (1956)-Finite Difference 
© Finite Element (15 iterations) of reasonable care and judgment in 
the apportionment and placement of 
, elements, even with seemingly 
First Trial Position of Free rather coarse configurations. Al- 
— though no comprehensive rules can 
———— be established with regard to the 
_— size and number of elements used 
r ne 1 in any given case, satisfactory con- 
(c) Comparison of Free Surface Positions as Caculated by Various Methods vergence was obtained for the prob- 
lems considered here by idealiza- 
Figure 3. Flow through a dam with vertical faces. tions ranging from 125 nodes and 
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TABLE 1 


COMPARISON OF FLOW PARAMETERS FOR SEEPAGE THROUGH A 
RECTANGULAR DAM WITH VERTICAL FACES 


Maximum H,/Hi 


Hi/He L Number of Variation in Q/kH, 
Iterations Free Surface Polubarinova- : 
Pressure FEM Kochina BENE Dupe 
0.100 20.0 15 +7 x.107° 0.1680 0.1143 0.2475 0.2475 
0.125 19.6 15 +6 x 107° 0.2821 - 0.2511 0.2511 
0.300 20.0 14 +3 x 1077 0.0002 0.0003 0.2275 0.2275 
0.500 20.0 8 +7 x 107° 0 0 0.1875 0.1875 
0.700 20.0 6 £7 x 107° 0 0 0.1275 0.1275 
0.900 20.0 5 +2 x 107° 0 ) 0.0475 0.0475 
Note: Hz = 10.0;k = 1.0; maximum iterations = 15; beta = 0.9; and tolerance = 0.001. 


101 elements to 461 nodes and 447 elements, and the number of iterations required to 
locate the free surface varied from 3 to 10. 


FLOW TO AN EXCAVATION IN A LAYERED AQUIFER SYSTEM 


For the case of flow to the type of excavation shown in Figure la, the following 
boundary conditions were used with the finite-element idealization. Zero normal flow 
conditions were imposed along the line of symmetry, the phreatic surface, and the im- 
permeable lower boundary; the boundary nodes at distance L were subjected to a hydro- 
static pressure distribution; zero pressure was specified at the nodes that describe 
the lower boundary of the excavation; and flow was allowed to occur through the excava- 
tion side wall. Finite-element meshes ranged in size from 125 nodes and 101 elements 
for the narrow trench with deep penetration to 236 nodes and 215 elements for the wide 
trench with shallow penetration. Relatively small, nearly square elements (2.5:2 units) 
were used in areas of rapidly changing pressures, whereas larger rectangular ele- 
ments (60:30 units) were used in areas with small pressure changes. Various sizes of 
rectangular elements were used in the transition regions. Because of the complex 
free-surface configurations that resulted from the various layered systems and perme- 
ability ratios, 15 different meshes were used for this problem, 

The number of variables needed to characterize the excavation problem shown in 
Figure la was reduced as follows: (a) the layer thicknesses D, were assigned a value 
of H,/3 equal to 30 units; (b) the effective length L of the domain was taken to be 1,000 
units; (c) the depth of penetration of the trench (H, - H,)/H, was chosen to be H,/6, 
H,/2, and 5H,/6; (d) the width W of the excavation was given values of 80 units and 
4 units to obtain information on width effects; (e) the excavation was assumed to be de- 
watered and to have permeable side walls; and (f) the relative permeability values 
(ratio of variable permeability k, to reference permeability k,) were taken to be 0.01, 
0.1, 1,10, and 100. Each layer, in turn, was considered to be the variable layer 
(whereas the permeability of the other 2 layers was held constant), and its relative 
permeability was allowed to traverse the full range of assumed values, thus yielding a 
2- or 3-layer system composed of 2 materials. This was done for each combination of 
excavation width and penetration, and information regarding the exit point on the 
seepage face and the flow quantity is given in Table 2. 


Interpretation of Results 


The location of the exit point varies as a function of the confined nature of the sys- 
tem, the exit point being higher when the uppermost layer (or layers) is less permeable 
than the underlying layer (or layers). This effect is present, to a degree, whether or 
not the excavation extends into the more pervious underlying layer, and it may also be 
observed when the layer that contains the bottom of the excavation is least permeable 
and the overlying layers are more permeable. In addition, for these same situations, 
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TABLE 2 
SUMMARY OF SEEPAGE CHARACTERISTICS FOR AN EXCAVATION 


Exit Point of 


Phreatic Surface Elon Parameter 


Relative Permeability 


Penetration Ratio of Layer 


4-Unit g0-Unit 4-Unit 80-Unit 


K/k,  ky/ky k3/Kp ~Wiath width. «Width «= Width 


% 1.00 1.00 1.00 0.83 0.83 0.169 0.184 
1.00 1.00 0.01 0.90 0.87 0.016 0.046 
1.00 1.00 0.10 0.85 0.83 0.076 0.116 
1.00 1.00 10.00 0.83 0.83 0.643 0.654 
1.00 1.00 100.00 0.83 0.83 5.288 5.336 
1.00 0.01 1.00 0.83 0.83 0.089 0.094 
1.00 0.10 1.00 0.83 0.83 0.111 0.118 
1.00 10.00 1.00 0.84 0.83 0.541 0.724 
1.00 100.00 1.00 0.89 0.85 1.436 3.450 
0.01 1.00 1.00 0.83 0.83 0.115 0.120 
0.10 1.00 1.00 0.83 0.83 0.120 0.126 

10.00 1.00 1.00 0.84 0.83 0.475 0.647 
100.00 1.00 1.00 0.86 0.84 0.995 2.080 

“a 1.00 1.00 1.00 0.50 0.50 0.149 0.156 
1.00 1.00 0.01 0,50 0.50 0.126 0.131 
1.00 1.00 0.10 0.50 0.50 0.128 0.133 
1.00 1.00 10.00 0.51 0.50 0.314 0.341 
1.00 1.00 100.00 0.89 0.85 1.404 1.721 
1.00 0.01 1.00 0.93 0.83 0.019 0.045 
1.00 0.10 1.00 0.53 0.50 0.071 0.089 
1.00 10.00 1.00 0.50 0.50 0.720 0.728 
1.00 100.00 1.00 0,50 0.50 6.416 6.466 
0.01 1.00 1.00 0,50 0.50 0.086 0.086 
0.10 1.00 1.00 0.50 0.50 0.092 0.093 

10.00 1.00 1.00 0.58 0.51 0.530 0.702 
100.00 1.00 1.00 0.77 0.63 1.534 3.327 

%s 1.00 1.00 1.00 0.17 0.17 0.119 0.119 
1.00 1.00 0.01 0,17 0.17 0.105 0.105 
1.00 1.00 0.10 0.17 0.17 0.106 0.107 
1.00 1.00 10.00 0.17 0.17 0.196 0.198 
1.00 1.00 100.00 0.86 0.85 0.996 1.071 
1.00 0.01 1.00 0.17 0.17 0.077 0.078 
1.00 0.10 1.00 0.17 0.17 0.082 0.083 
1,00 10.00 1.00 0.31 0.30 0.442 0.445 
1,00 100.00 1.00 0.84 0.83 1.372 1.517 
0.01 1,00 1.00 0.91 0.90 0.015 0.017 
0.10 1.00 1.00 0.45 0.42 0.060 0.062 

10.00 1.00 1.00 0.17 0.17 0.714 0.714 
100.00 1.00 1.00 0.17 0.17 6.650 6.654 


the case with the narrower excavation width indicates consistently higher values for 
the exit point elevation. 

For a penetration ratio, (H, - H,)/H,, of 1/6, the results of which are shown in 
Figures 4a and 4b, the relative permeability of layer 3 exerts the most influence on 
the flow quantity, and little variation in flow quantity is caused by changes in k, andk, 
when their relative permeability values are approximately less than 1. In addition, the 
effect of the excavation width can be evaluated by comparing the 2 sets of curves. In 
the case of the wide excavation, there is very little variation in flow quantity with a 
variation in any of the layer relative permeabilities, except for the very high or very 
low relative permeability values; whereas, in the case of the narrow excavation, the 
influence of layer 3 manifests itself more readily. Another quantitative appraisal of 
excavation width can be obtained by comparing the flow quantities for a given set of 
permeability conditions; the difference between any two such values essentially repre- 
sents the additional quantity of flow that is passing through the bottom of the excava- 
tion. The curves for layer 1 indicate that it exerts the least influence on the flow char- 
acteristics of the system. 

As the depth of the excavation extends into layer 2, that layer exerts the largest 
influence on the flow in the system, and these effects are shown in Figures 4c and 4d 
for the different excavation widths. Figure 4c shows that the underlying layer (layer 1) 
has a dominant influence in the relative permeability range from 1 to 10, even though 
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it has not been penetrated. Also, a reduction in the influence of layers 1 and 3 may be 
observed for relative permeabilities of more than 1, and an increase in the influence of 
layer 2 for relative permeabilities of less than 1 is apparent; these changes are attrib- 
utable to excavation width. 

Figures 4e and 4f show the results for the case where the excavation penetrates the 
bottom layer. For this situation, layer 1 exerts the major effect on the flow quantity. 
A comparison of these results shows very little change in flow due to excavation width 
because the underlying layer is impermeable. 


Conclusions 


Several qualitative conclusions may be drawn from the results of flow studies. For 
the geometrical configuration studied, it is apparent that the layer containing the bot- 
tom of the excavation manifests the dominant influence on the flow quantity. Where the 
excavation is wide and an underlying layer is relatively permeable, the underlying layer 
may, for a certain range of relative permeabilities, contribute more to the flow quantity 
than the layer that contains the bottom of the excavation. The exit point of the phreatic 
surface may be expected to assume a location that is a function of the degree to which 
the system is confined. In particular, there are certain commonly encountered field 
conditions that make it difficult to lower the free surface. 


SUBDRAINS IN LAYERED AQUIFERS 


In a more specific application of the preceding concept, the cross sections shown in 
Figure 1b were considered (18), and the effect of varying the thickness of layer 2 on the 
flow characteristics for several subdrain penetration depths was investigated. Although 
this case represents a subdrain system that is quite general in nature, the particular 
combination of soil profile, varying layer thicknesses, and relative permeabilities is 
felt to typify a portion of the proposed Crosstown Expressway in Chicago. Before the 
specific results obtained in this phase of the study are considered, several deductions 
may be made from the preceding results. First of all, for the particular combination 
of relative permeabilities shown, it should prove difficult to lower the phreatic surface 
sufficiently unless the silt-sand layer is tapped directly by the drains. Second, changes 
in the spacing of the subdrains will probably have little effect on the resulting flow 
quantities (analogous to the width effect for the excavation) because most practical 
spacings for highway drains will approximate a wide excavation. Third, changes in the 
thickness of the silt layer will probably not have a significant effect on the quantity of 
flow until the subdrains tap the silt layer or until it becomes very thick (i.e., the sand 
layer will maintain an almost constant pressure on the lower interface of the silt 
layer). 

With the foregoing thoughts in mind, we obtained solutions for drain spacing widths 
W of 50 units and 200 units and for penetration ratios of 1/6, 1/3, 1/2, and 1, where the 
water level in the drains was assumed to be at the penetration depth. Portions of the 
resulting data are given in Table 3 and shown in Figure 5. Table 3 gives information 
that indicates the position of the phreatic surface (offset values from the original water 
table are shown for various points along the effective length L of the domain). The 
value for L was 10 times the penetration ratio (H, - H,), each drain was 2 units wide, 
the thickness D, of the silt layer was set at 0, 1/3, 2/3, and 1 times D, whereas D, was 
a constant 30 units, and H, was held constant at 60 units. 

For the purpose of this study, the absolute magnitudes of the flow quantities were 
considered to be of little concern. The major interest centered around the effect of the 
thickness of layer 2 on the relative flow quantities and the location of the free surface. 
Thus, the case where D, equals zero was selected as a reference case for each differ- 
ent drain penetration, and the flow quantity Q for various values of D, is divided by the 
flow quantity Q, to form a dimensionless flow parameter Q/Q,, which is shown in Fig- 
ure 5 as a function of the thickness ratio D,/D, and the depth of penetration parameter 
(H, - H,)/H, When D,/D, is zero, the bottom layer is all sand, and the corresponding 
value of Q/Q, is 1 for all cases. When D,/D, equals 1, the bottom layer is all silt, and 
appropriate values for Q/Q, are shown in Figure 5. 


TABLE 3 
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LOCATION OF PHREATIC SURFACE FOR FLOW TO DRAINAGE TRENCHES 


50-Unit Width 


200-Unit Width 


P i Layer 
enetration Thickness 
Ratio Ratio Line of 1Unit 1Unit 0.005L 0.40L Line of 1Unit i1Unit 0.05L 040L 
Symmetry Left Right Right Right Symmetry Left Right Right Right 
Ye 1 5.25 9.80 9.00 7.10 2.27 3.16 8.60 8.70 7.12 2.30 
7h 3.83 8.05 8.05 6.24 0.96 0.12 8.05 8.20 6.22 0.95 
hs Beal 8.00 8.00 6.18 0.89 0.09 8.00 8.10 6.16 0.88 
0 3.58 8.00 8.00 6.11 0.83 0.07 8.00 8.10 6.11 0.83 
a 1 12.89 17.70 17.20 13.37 6.67 11,27 17.50 17.40 13.41 6.66 
* 7.13 13.10 11.60 7.80 0.29 0.26 11.30 11.60 7.79 0.29 
hs 6.75 12.80 11.50 7.53 0.22 0.17 11.00 11.20 7.54 0.22 
0 6.38 12.10 11.10 7.28 0.18 0,13 10.80 10.80 7.28 0.18 
%a a 28.44 29.90 29.60 26.97 16.10 28.28 29.80 29.60 26.98 16.51 
WA 12.77 18.50 17.00 9.35 0.72 1,49 12.00 14.00 9.37 0.34 
/s 9.05 13.80 13.70 7.02 0.49 0.94 13.50 12.80 7.04 0.49 
0 28.47 30.00 29.50 27.06 16.83 28,31 30.00 30.00 27.07 16.83 
1 1 60.00 60.00 57.30 48.22 26.51 60,00 60.00 57.30 48.22 26.51 
7h 60.00 60.00 57.40 50.04 30.00 60.00 60.00 57.40 50.04 30.00 
hs 60.00 60.00 57.60 52.51 34,00 60.00 60.00 57.60 52.51 34.00 
0 60.00 60.00 57.10 48.43 27,31 60.00 60.00 57.10 48.43 27,31 


The data given in Table 3 show that the spacing of the drains has little effect on the 
position of that portion of the free surface located outside the drains. However, the ex- 
tent of the drawdown between the drains is significantly affected when the underlying 
layers are a combination of sand and silt but not when only 1 material is present for a 
penetration ratio of 1/2. The effect of drain spacing on flow quantity is negligible for 
all cases, and the curves shown in Figure 5 represent both values of W. Although the 
thickness of the silt layer has relatively little effect on flow quantity when the drains 
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Figure 5. Relative flow quantity versus layer thick- 
ness ratio for various drain penetrations. 


do not tap that layer, a large effect is ob- 
served when the drains do tap the silt layer. 
However, full penetration of the drainage 
trenches reduces the effect of the silt layer, 
because the sand aquifer contributes the 
major portion of the flow quantity for all 
cases, 


Results of Subdrain Studies 


Although the spacing of drains in an open 
cut has been shown to have little effect on 
flow quantities and on the position of that 
portion of the free surface located outside 
of the drains, it does exert a substantial in- 
fluence on the drawdown between the drains. 
However, under certain conditions the thick- 
ness of an intermediate silt layer, typical of 
a situation encountered along a portion of a 
proposed highway cut in the Chicago area, 
was found to affect significantly the seepage 
characteristics of the system. The results 
obtained from this latter, more specific 
study are qualitatively consistent with those 
that were deduced from the preceding, more 
general study of flow into an excavation. 
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APPLICATION OF QUALITATIVE RESULTS TO AN EXAMPLE PROBLEM 


The area to be traversed by the proposed Crosstown Expressway in Chicago consists 
of various surficial soils, lacustrine sands and silts, glacial till sheets, pockets of 
dense granular soils, and bedrock. Conditions are quite variable, and the water table 
is approximately 5 to 10 ft beneath the existing ground surface in most places. Ac- 
cordingly, depressed sections of the roadway will generally be located in soil profiles 
very similar to those considered in the previous sections of this study. As a particu- 
lar example, consider the typical soil profile (obtained for the preliminary subdrain 
study) shown in Figure 6a with its representative permeability values. Because a de- 
pressed roadway is to be constructed in this profile, it is necessary to specify the type 
and placement of permanent drains to maintain satisfactory dewatering and free-surface 
location. 

Based on the preceding study, the system appears to display "confined flow" char- 
acteristics, and hence the free surface will be difficult to lower. Narrow drainage 
trenches will probably be used, and therefore no increase in flow quantity should be 
expected for each drain because of width effect. However, the use of multiple trenches 
will create a pseudo-wide-trench response, and disproportional amounts of water will 
be drawn from the underlying strata by creating a large region in which vertical flow 
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Figure 6. Typical soil profile and drainage systems for a roadway cut in the 
Chicago Crosstown Expressway. 
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predominates. Furthermore, it would appear that at least 2 drains are necessary, 
one on either side of the roadway, and that they must penetrate the layer that has a 
10-5 cm/sec permeability to lower the free surface. Also, these drains should be 
spaced as close together as possible to maximize the drawdown between the drains. 

The preceding qualitative deductions were checked by studying 4 different drain 
configurations by means of the finite-element method. Idealizations involving from 
236 nodes and 223 elements to 265 nodes and 252 elements were used, and between 3 
and 7 iterations were needed to obtain the free-surface location. Although Figure 6b 
shows what might be a first intuitive selection of drain penetration, the results of the 
earlier studies presented here have demonstrated that a configuration of this type 
cannot provide satisfactory drainage. This conclusion was verified by the finite- 
element solution, which is shown in Figure 6b. As can readily be seen, the phreatic 
surface is totally unsatisfactory; and problems with slope stability, subbase and retain- 
ing wall drainage, and excess pressures on the pavement will be encountered if no 
supplementary drains are supplied. 

Extension of the drains into the more pervious silt layer, the next logical choice in 
light of the general study, yields the satisfactory solution shown in Figure 6c, as long 
as a minimum spacing of drains is maintained. According to the preceding study, just 
tapping the silt layer will give essentially the same position for that portion of the free 
surface outside the drains and the same flow quantity into the drains. However, the 
position of the free surface between the drains will need to be controlled by use of a 
third or centerline drain, as shown in Figure 6d. This latter solution, which was veri- 
fied by a finite-element computation, has the advantage of requiring smaller quantities 
of excavation and filter materials. 

In each of these cases, no special provision for drainage of the retaining wall was 
considered. The calculated seepage quantity for the final drainage system shown was 
less than 2 cu ft per day per foot of length for each drain, with the bottom layer con- 
tributing approximately 96 percent of this quantity. Thus, the qualitative predictions, 
based on the results of previous studies, for the seepage characteristics of a typical 
depressed roadway in a layered system provided a reasonably accurate description of 
the results that were obtained from a finite-element analysis. 


SUMMARY 


The qualitative concepts developed for the general case of an excavation ina layered 
aquifer are shown to apply to more specific field problems that involve the installation 
of subdrains. These data aid in understanding the complex interaction that occurs in 
plane flow through systems of layered aquifers, and they allow a more direct and 
logical approach to the design of highway drainage systems. 
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